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Where a translation of the intemational application must be furnished to an elected Office, that translation must 
contain a translation of any annexes to the intemational preliminary examination report. It is the applicant's 
responsibility to pTSpare and furnish such translation directly to each elected Office concemed. 

For further details on the applicable time limits and r«qulramsnts of the alectad Offices, see Volume II of the 
PCT Applicant's Guide. 



Namo end mailing address ot mo IPEA/ 

.™ European Patent Ottk» 

mii 0^298 Munich 

Tel, +49 89 2399 - 0 T^: S238S6 Opmu d 
Pax: +49 89 2399 - 446S 



Authorized officer 
Schmethusen. S 

Tel. +49 89 2399^15*?/ 7^^? >f 




Form PCT/lPeA/416 (July 1992) 



30-OB-QQ 12:20 From: HARRISON GODDARD FOOTE 



+441142730312 



T-663 P. 99/99 Job-687 



PATENT COOPERATION TREATY 

PCT 

INTERNATIONAL PRELIMINARY EXAMINATION REPORT 

(PCT Article 36 and Rule 70) 



Applicant's or a0ant!& (ila referenca 


Soe Notification ai Transmittal ot International 


P30S24WO 


FOR FURTHER ACTION Praliminary Examination Report (Fofm PCT/lP£A;4t6) 


tntamational applicatkin No. 


Iniarnatjonai filing date (day/monih/year) 


Priority date (day/month/y9ar) 


PCT/GB99/00550 


02/03/1999 


03/03/1998 


Intsfnational Patent aassification (IPC) or national classification and IPC 




H03H 17/02 






Applicant 






UNIVERSITY OF SHEFFIELD et a|. 







1 . This international preliminary examination report has been prepared by this International Preliminary Examining Authority 
and is transmitted to the applicant according to Article 36, 

2. This REPORT consists of a tcrtal of 5 ahoots, including this cover sheet. 

® This report is also accompanied by ANNEXES, i.e. sheets of the description, claims and/or drawings which have 
been amended and are the basis for this report and/or sheets containing rectifications mads before this Authority 
(see Role 70.16 and Section 807 of the Administrative Instructions under the PCT). 

These annexes consist of a total of 3 sheets. 



3. This report contains indications relating to the following items; 
I S Basis of the report 



II 
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Priority 


III 
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Non-establlshment of opinion with regard to novelty, inventive step and industrial applicability 


IV 




Lack of unity of invention 


V 




Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citetions and explanations siiporting such statement 


VI 




Certain^cumants cited 


VII 
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Certain dafacta in the internationat application 


VIII 
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Certain observations on the international application 
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EXAMINATION REPORT 



International application No. PCT/GB99/00550 



I. Basis of the report 

1 This report has been drawn on the basis of (substitute sheets which have been furnished to the receiving Office in 
response to an invitation under Article 14 are referred to in this report as 'originally filed" and are not annexed to 
the report since they do not contain amendn)ents.): 

Description, pages: 

1-78' as originally tiled 



Claiina, No.: 

1 *26.27 (pan) as originally filed 

27 (part), 28-32 as receh/sd on 



07/03/2000 with letter of 



07/03/2000 



Drawings, sheets: 

1 /38-38/3a as orioinaliy filed 



2. The amendments have resulted in the cancellation of: 

□ the description, pages; 

a the claims, Nos,: 33^7 

n the drawings, sheets; 



3. □ This report has been established as it (soma of) the amendments had not been made, since they have been 
considered to go beyond the disclosure as filed (Rule 70,2(c)): 



4- Additional obsen/ations, it necessary: 
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EXAMINATION REPORT International application No. PCT/GB99/QQ55Q 

V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial 
applicability; citations and explanations supporting such statement 

1. Statement 
Novelty (N) 

Inventive step (IS) 



Yes: 


Claims 


1-29.31,32 


No: 


Claims 


30 


Yes: 


Claims 


1-29,31,32 


No: 


Claims 


30 


Yes: 


Claims 


1-32 


No: 


Claims 





2. Citations and explanations- 
see separate sheet 

VI. Certain documents cited 

1. Certain published documents (Rule 70.10) 
and / or 

2. Non-written disclosures (Rule 70,9) 
see separate sheet 
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international application No. PCT/GB99/00550 



EXAMINATION REPORT - SEPARATE SHEET ^ ^ 

TO SECTION V 

1 . The first 3 documents cited in the search report are not considered particularly 
relevant to the subject matter disclosed in the present application. However the 
paper "Spectral Amplitude Warping (SAW) for noise spectrum shaping in audio 
coding". R. Lefebvre, Laflamme, 1997 IEEE ICASSP. April 1997, pp 335-338 is 
considered to be the closest prior art. This document discloses a non linear 
transformation applied to each frequency sample of the discrete spectrum of a 
speech signal. 

2. Independent claims 

2.1 The subject matter of claim 1 (design method) and 1 1 (equivalent design system), 
an in particular the transfer of energy between two specified spectra using the 
coefficients of a time or spatial domain description of a generalised non-linear 
system, is neither disclosed nor rendered obvious by the available state of the art 
and is therefore considered to fulfil the requirements of Article 33(2) and (3) PCT. 

The same applies to claim 10, which contains, by reference, the subject matter of 
claim 1- 

The same applies also to claim 20, claiming a computer program product for 
realizing steps identical to those of claim 1 , and to claim 32 claiming an article of 
manufacture comprising such a computer program. 

The sarTie applies also to claim 29, claiming a non-linear system based on the 
design metltod of claim 1 . 

The same applies to claim 31 . comprising a data processing system as claimed in 
any of claims 11 to 19. 

2.2 The subject matter of claim 30 recites a non-linear system for transferring energy, 
which must is interpreted as meaning merely "suitable for transferring energy". 
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and characterized by comprising means for implementing a method as claimed in 
any of claims 1 to 9, again interpreted as meaning merely "means suitable for 
implementing" such methods. Therefore, claim 30 is not considered novel 
because it does not recite any means which unambiguously and specifically 
contribute towards executing the novel and inventive method of claim 1 , and in its 
present form amounts to a simple computer. 

3. Dependent claims 

The dependent cJaims 2-9,12-19,21-28 are considered novel and inventive for the 
same reason as the independent claim they depend from. 



TO SECTION VI 

1 . The paper "Narx Approach to black-box modelling of circuit elements", by I. A 
Maio, I.S. Stievano, F.G. Canavero, Proceedings of the 1998 IEEE International 
Symposium on Circuits and Systems, pp 41 1-414, Monterey, 31 May-03 June 
1998 was made available to the public between the priority date of the present 
application (03 March 1998) ane its filing date (02 March 1999). 
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the description for the system output spectrum given in 
terms of the spectrum of the input signal and the 
coefficients of the time or spatial domain model of the 
said generalised non^linear system. 



28, 



10 



15 



20 



25 



30 



A computer, program product as claimed in any of claims 
20 to 21, further comprising computer program code 
means for substituting H„( jw,,--^, jwj given in (C21) into 
<C20), and substituting the resultant expression for 
Y„(jw) into (C19) 'to obtain the description for the 
system output spectriom in terms of a function of the 
■spectrum of the input signal and the coefficients of 
the- time or spatial domain model of the generalised 
nonlinear system. 



29, A method for designing a non-linear system for 
^transferring energy from a time or spatial domain input 
signal having a first spectrum at a first pre- 
determinaiDle frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinabie frequency or 
range of frequencies substantially as described herein 
with reference to and/or as illustrated in the 
accompanying drawings . 



30. A method for manufacturing a non-linear system 
comprising the steps of designing a non-linear system 
for transferring energy from a time or spatial domain 
input signal having a first spectrum at a first pre- 
determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies using a method substantially as 
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described herein with reference to and/or as 
illustrated in the accompanying drawings and materially 
producing the non-linear system so designed. 



5 31. A data processing system for designing a non-linear 
system 

for transferring energy from a time or spatial domain 
input signal having a first spectrum at a first pre- 
determinable frequency or range of frequencies to a 
10 time or sparial domain output sigrial having a second 

spectrum at a second pre-determlnable frequency or 
■ -range of- frequencies substantially as described herein 
with reference to and/or as illustrated in the 
accompanying drawings . 

15 

32. A computer program product comprising computer program 
'code means for designing and/or implementing a non- 
linear system for transferring energy from a time or 
spatial domain input signal having a first spectrum at 

20 a first prd-determinable frequency or - range of 

frequencies to a time or spatial domain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies substantially as 
described herein with reference to and/or as 

25 illustrated in the accompanying drawings , 

33. A non-linear system for transferring energy from a time 
or spatial domain input signal having a first spectrum 
at a first predeterminable frequency or range of 
frequencies to a time or spatial domain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies substantially as 
described herein with reference to and/or as 
illustrated in the accompanying drawings . 



30 
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34. 



10 



13 



25 



30 



35. 



A non-linear system which can transfer energy from a 
time or spatial domain input signal having a first 
spectrum at a first pre-determinabie frequency or range 
of frequencies to a time or spatial domain output 
signal having a second spectrum at a second pre- 
determinable frequency or range of f re<5uencies , said 
system comprising 

means for identifying the first spectrtim of the time or 
spatial domain inpup signal from which energy is to be 
transferred, 

means for specifying the second spectrum of the time or 
spatial domain output signal to which said energy is to 
be transferred, and 



means for giving effect to the energy transfer using 
coefficients of a time or spatial domain description of 
a generalised non-linear system, said coefficients 
having been calculated using a frequency domain 
20 description of said output signal, for example, the 

output spectrum, expressed in terms of a frequency 
domain description of said input signal and 
coefficients of a time or spatial domain description of 
a generalised non-linear system. 



A non-linear system for transferring energy from a time 
or spatial domain input signal having a first spectrum 
at a first pre-determinabie frequency or range of 
frecjuencies to a time or spatial domain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies, said system 
comprising means for implementing a method as alaimed 
in any of claims 1 to- 9 or 29. 



35 
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36. A non-linear system for transferring, energy from a time 
or spatial domain input signal having a first spectrum 
at a first pre-determinable frequency or range of 
frequencies to a time or spatial domain output signal 
having a second spectrum it a second pre-determinable 
frequency or range of frequencies, said noa-lineax 
system comprising a data processing system as claimed 
in any of claims 11 to 19 or 31. 



An article of manufacture comprising a computer usable 
medium with computer readable program code means 
einbodied in che medium for designing a non-linear 
3.ystem-for transferring energy from a time or spatial 
domain input signal having a first spectrum a first 
pte-determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies, the computer readable program 
coda means in said article comprising: 

computer readable prograin code means for identifying, 
the first apectruin of a time or spatial domain input 
signal from which energy is to be transferred, 
computer readable prograin code means for specifying the 
second spectrum of a time or spacial domain output 
signal to which said energy is to be transferred, and 
<;emputer readable program code means for calculating, 
usin« a frequency domain description of the output 
signal, for example, the output spectrum, expressed in 
terms of a frequency domain description of the input 
30 and coefficients of a time or spatial domain 

description of a generalised non-linear system, the 
coefficients of a time or spatial domain description of 
said generalised non-linear system in order to give 
effect to the energy transfer. 
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Applicant's or agent's file reference 


FOR FURTHER ®®® Notification of Transmittal of International Search Report 

(Form PCT/ISA/220) as well as, where applicable, item 5 below. 

ACTION 


International application No. 

PCT/GB 99/00550 


International filing date (day/month/year) 

02/03/1999 


(Earliest) Priority Date (day/month/year) 

03/03/1998 


Applicant 

UNIVERSITY OF SHEFFIELD et al . 



This International Search Report has been prepared by this International Searching Authority and is transmitted to the applicant 
according to Article 18. A copy is being transmitted to the International Bureau. 

This International Search Report consists of a total of 3 sheets. 

[X] It is also accompanied by a copy of each prior art document cited in this report. 



1 . Basis of the report 

a. 



With regard to the language, the international search was carried out on the basis of the international application in the 
language in which it was filed, unless otherwise indicated under this item. 



□ 



the international search was carried out on the basis of a translation of the international application furnished to this 
Authority (Rule 23.1(b)). 



With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international search 

was carried out on the basis of the sequence listing : 

I I contained in the international application in written form. 

filed together with the international application in computer readable form. 

furnished subsequently to this Authority in written form. 

furnished subsequently to this Authority in computer readble form. 



□ 
□ 
□ 
□ 



V 



2. 
3. 



□ 

□ 
□ 



the statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 
international application as filed has been furnished. 

the statement that the information recorded in computer readable form is identical to the written sequence listing has been 
furnished 

Certain claims were found unsearchable (See Box I). 
Unity of invention is lacking (see Box II). 



4. With regard to the title, 

[X| the text is approved as submitted by the applicant. 

I I the text has been established by this Authority to read as follows: 



5, With regard to the abstract, 

PC] the text is approved as submitted by the applicant. 

□ the text has been established, according to Rule 38.2(b), by this Authority as it appears in Box III. The applicant may, 
within one month from the date of mailing of this international search report, submit comments to this Authority. 

6. The figure of the drawings to be published with the abstract is Figure No. 4 

I I as suggested by the applicant. Q None of the figures. 

pC] because the applicant failed to suggest a figure. 

I I because this figure better characterizes the invention. 
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"O" document referring to an oral disclosure, use, exhibition or 
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"P" document published prior to the international filing date but 
later than the priority date claimed 



'T" later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underiying the 
invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
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I. Basis of the report 

1 . This report has been drawn on the basis of {substitute sheets which have been furnished to the receiving Office in 
response to an invitation under Articie 14 are referred to in this report as "originaiiy filed" and are not annexed to 
the report since they do not contain amendments,): 

Description, pages: 

1 -78' as originally filed 

Claims, No.: 

1 -26,27 (part) as originally filed 

27 (part),28-32 as received on 07/03/2000 with letter of 07/03/2000 

Drawings, sheets: 

1/38-38/38 as originally filed 

2. The amendments have resulted in the cancellation of: 

□ the description, pages: 

K the claims, Nos.: 33-37 

□ the drawings, sheets: 

3. □ This report has been established as if (some of) the amendments had not been made, since they have been 

considered to go beyond the disclosure as filed (Rule 70.2(c)): 

4. Additional observations, if necessary: 
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V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial 
applicability; citations and explanations supporting such statement 

1 , Statement 



Novelty (N) 


Yes: 


Claims 


1-29,31,32 




No: 


Claims 


30 


Inventive step (IS) 


Yes: 


Claims 


1-29,31,32 




No: 


Claims 


30 


Industrial applicability (lA) 


Yes: 


Claims 


1-32 




No: 


Claims 





2. Citations and explanations 
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1. Certain published documents (Rule 70.10) 
and / or 

2. Non-written disclosures (Rule 70.9) 
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TO SECTION V 



1 . The first 3 documents cited in the search report are not considered particularly 
relevant to the subject matter disclosed in the present application. However the 
paper "Spectral Amplitude Warping (SAW) for noise spectrum shaping in audio 
coding", R. Lefebvre, C. Laflamme, 1997 IEEE ICASSP, April 1997, pp 335-338 is 
considered to be the closest prior art. This document discloses a non linear 
transformation applied to each frequency sample of the discrete spectrum of a 
speech signal. 

2. Independent claims 

2.1 The subject matter of claim 1 (design method) and 1 1 (equivalent design system), 
an in particular the transfer of energy between two specified spectra using the 
coefficients of a time or spatial domain description of a generalised non-linear 
system, is neither disclosed nor rendered obvious by the available state of the art 
and is therefore considered to fulfil the requirements of Article 33(2) and (3) PCT. 

The same applies to claim 10, which contains, by reference, the subject matter of 
claim 1. 

The same applies also to claim 20, claiming a computer program product for 
realizing steps identical to those of claim 1, and to claim 32 claiming an article of 
manufacture comprising such a computer program. 

The same applies also to claim 29, claiming a non-linear system based on the 
design method of claim 1 . 

The same applies to claim 31, comprising a data processing system as claimed in 
any of claims 11 to 19. 

2.2 The subject matter of claim 30 recites a non-linear system for transferring energy, 
which must is interpreted as meaning merely "suitable for transferring energy", 
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and characterized by comprising means for implementing a method as claimed in 
any of claims 1 to 9, again interpreted as meaning merely "means suitable for 
implementing" such methods. Therefore, claim 30 is not considered novel 
because it does not recite any means which unambiguously and specifically 
contribute towards executing the novel and inventive method of claim 1, and in its 
present form amounts to a simple computer. 

3. Dependent claims 

The dependent claims 2-9,12-19,21-28 are considered novel and inventive for the 
same reason as the independent claim they depend from. 

TO SECTION VI 

1. The paper "Narx Approach to black-box modelling of circuit elements", by LA 
Maio, I.S. Stievano, F.G. Canavero, Proceedings of the 1998 IEEE International 
Symposium on Circuits and Systems, pp 411-414, Monterey, 31 May-03 June 
1998 was made available to the public between the priority date of the present 
application (03 March 1998) ane its filing date (02 March 1999). 
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takir.g N as N and determining the coefficients of the 
time or spatial dorciaiu model cf che generalised non- 
linear system from the descriprion for the system 
output spectr^um given in terrt^s of the spectrum of the 
5 input signal and the coefficients of th© time or 

soatxal domain model of the said generalised non-linear 
system- 



28- A computer program product as claimed in any of claiias 
10 20 to 21 r further comprising computer program code 

means for substituting jw^,-", j v/^) given in s'C2X) into 
(C20), and substituting the resultant expression for 
Y^(jw) into {C19) to obtain the description for the 
c^vfT^-'-m output spectrum in ne-rms. of a functio- o^^ rhe 
15 spectrum of the input signal and the coefficients of 

the time or spatial domain m.odel of the generalised 
nonlinear system. 
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25 



30 



29. A non-linear system which can transfer energy from a 
time or spatial domain input signal having a first 
spectrum at a first pre-determinable frequency or range 
of frequencies to a time or spatial domain output 
signal having a second spectrum at a second pre- 
determinatile frequency or range of frequencies, said 
system comprising 

means for identifying the first apectr\jrti of tHe time or 
spatial domain input signal from which energy is to be 
transferred, 

means for specifying the second spectrum of the time or 
spatial domain output signal to which said energy is to 
be transferred, and 
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25 



means for giving effect to rhe energy traR.sfer usir.g 
corfficients of a time or spacial domain description, of 
a generalised non-linear system, said coefficients 
havia'3 been calculated using a frequency domain 
deacriptiorv o£ said outputi signal, for example, th.© 
output spectruir., expressed in terrrs of a frequency 
docTvain description of said input signal and 
coefficients of a time or spatial domain description cf 
a generalised non-linear system. 



30. Px non-Linear systsm to transferring energy from a time 
or spatial domain input signal having a first spectrum 
at a first pre-determinable frequency or range of 
'requencies cc vine cr sp^"ial dot^ain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies, said system 
comprising means for implementing a method as claimed 
in any of claims 1 to 9 . 

31, A non-linear system for transferring energy from a rime 
or spatial domain input signal having a first spectrum 
at a first pre-determinable frequency or range of 
frequencies to a rime or spatial domain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies, said non-linear 
system comprising a data processing system as claimed 
in any of claims 11 to 19, 



30 32 . An article of manufacture comprising a computer usable 
medium with computer readable program code means 
embodied in the medium for designing a nan- linear 
system for transferring energy from a time or spatial 
domain input signal having a first spectrijm at a fxrst 

35 pre-determinable frequency or range of frequencies to a 
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time or spatial domain output sigr.ai having a second 
spectrura at a second pre-detsrini:i.abie frequency or 
range of f reqpaencies/ the computer readable program 
code means in said article comprising: 

cocnputer readable program code means for identifying 
the first spectrum of a time or spatial domain input 
signal from which energy is to be transferred/ 

computer readable program code means for specifying the 
second spectrum of a time or spatial domain outpuc 
signal to which said energy is to bs transferred^ and 

computer readable program code meens for c;a.lculat ing/ 
using a frequency domain description of the output 
signal, for example, uhe output spectrum, expressed in 
terms of a frequency domain descr-Lc-tior- of thr- : npi^-^ 
and coefficients of a time or spatial domain 
description of a generalised non-linear system, the 
coefficients of a time or spatial domain description of 
said generalised non-linear system in order to give 
effect to the energy transfer. 
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(57) Abstract 



The present invention relates to nonlinear systems and methods of design thereof in the frequency domain. Typically, conventional 
linear filter design involves attenuating signals at frequencies which are not of interest and dissipating the energy at those frequencies as, 
for example, heat or sound. However, in most sytems, it is not always convenient to design a linear system or design a system solely with 
energy attenuation in mind. Therefore, the present invention provides a nonlinear system and method of design thereof in the frequency 
domain which can be used to transfer energy at a first pre-determinable frequency, or frequency range, to a second pre-determinable 
frequency, or frequency range. Using the method of the present invention a nonlinear system can be developed which can meet given 
energy transfer requirements or a nonlinear system can be designed which can alter the transfer function of an existing nonlinear or linear 
system. 
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Nonlinear System^ Method o£ Design thereof and 
Con^duter Program Product 

The present invention relates to non-linear systems, 
5 methods of design thereof in the frequency domain and 
computer program products. More particularly, the present 
invention relates to a non-linear system having pre- 
deterrainable frequency response characteristics. The 
invention can be utilised to design and realise, for 
10 example, nonlinear filters having a required frequency 
response or transfer functions having specified transfer 
characteristics or within a control system context. 

The possible frequency components in an output signal 
15 of a linear system are exactly the same as the frequency 
components of a corresponding input signal. Conventional 
linear filter design is based on the principle that energy 
in unwanted frequency bands is attenuated. 



20 The Dolby filter, which varies the amplitude of the 

output signal as a function of the level and frequency of 
the input, is an example of a nonlinear filter system. 
However, when compared with the input, the output does not 
contain any additional frequency components. Modulation is 

25 another concept related to nonlinear filtering, which is 
associated with signal transmission where the signal to be 
transmitted is modulated by a carrier signal and then 
transmitted through a medium. Although a modulation device 
allows energy to be moved from one frequency band to 

30 another, the output frequency components of such a device 
depend not only on the input components but mainly on the 
carrier signal. Therefore, the energy transfer implemented 
by modulation is realised by a two input and one output 
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system where one input is the carrier signal and the other 
input is the signal to be processed. 

The prior art lacks a non-linear system and method 
5 of/apparatus for the design of such a non-linear system for 
predictably transferring energy from one frequency or band 
of frequencies of an input signal to another frequency or 
band of frequencies independently of any other input signal. 
Further, the prior art lacks a non-linear control system 
10 which can predictably transfer energy from one frequency 
band to another frequency band. 

It is an object of the present invention to at least 
mitigate some of the problems of the prior art. 

15 

Accordingly, a first aspect of the present invention 
provides a method for designing a non-linear system for 
transferring energy from a time or spatial domain input 
signal having a first spectrum at a first pre-determinable 
20 frequency or range of frequencies to a time or spatial 
domain output signal having a second spectriim at a second 
pre-determinable frequency or range of frequencies . 

Preferably, the method comprises the steps of 

identifying the first spectrum of the time or spatial domain 
25 input signal from which energy is to be transferred, 

specifying the second spectrum of the time or spatial domain 
output signal to which said energy is to be transferred, and 

calculating, using a frequency domain description of said 
output signal, for example, the output spectrum, expressed 
30 in terms of a frequency domain description of said input 
signal and coefficients of a time or spatial domain 
description of a generalised non-linear system, the 

2 
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coefficients of a said time or spatial domain description of 
said generalised non-linear system in order to give effect 
to the energy transfer. 



allows energy at a particular frequency within a given 
system to be transferred to another frequency or band of 
frequencies at which the response of the system is greatly 
reduced or negligible, or 

10 allows energy of a signal which is transmitted at a 
particular frequency or band ' of frequencies to be 
transferred, without using an additional modulating signal, 
to another frequency or band of frequencies at which the 
associated transmission media allows signals to pass, or 

15 allows energy at a particular band of frequencies to be 
transferred and spread over a new wider range of frequencies 
so as to attenuate the energy by employing the desired 
interkernel and intrakernel effects of nonlinear systems. 

20 The present invention is based upon the relationship 

between the input and output spectra or frequency components 
of nonlinear systems, and the relationship between the input 
and output frequencies and/or frequency ranges in the 
nonlinear case. In addition, the invention utilises a 

25 mapping between the time or spatial domains and frequency 
domain which allows the output spectra or frequency content 
of nonlinear systems to be described completely by the 
coefficients of time or spatial domain models which 
represent the filter or non-linear system to be constructed. 

30 A second aspect of the present invention provides a 

method for manufacturing a non-linear system for 
transferring energy from a time or spatial domain input 



5 



Advantageously, the present invention 



3 
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signal having a first spectrum at a first pre-determinable 
frequency or range of frequencies to a time or spatial 
domain output signal having a second spectrum at a second 
pre-determinable frequency or range of frequencies, said 
5 method of manufacture comprising the steps of 

(a) designing said non-linear system comprising the steps of 

identifying the first spectrum of the time or spatial domain 
input signal from which energy is to be transferred, 

specifying the second spectrum of the time or spatial domain 
10 output signal to which said energy is to be transferred, and 

calculating, using a frequency domain description of said 
output signal, for example, the output spectrum, expressed 
in terms of a frequency domain description of said input 
signal and coefficients of a time or spatial domain 
15 description of a generalised non-linear system, the 
coefficients of a time or spatial domain description of said 
generalised non-linear system in order to give effect to the 
energy transfer, and 

(b) materially producing the non-linear system so designed. 

20 

A third aspect of the present invention provides a data 
processing system which can transfer energy from a time or 
spatial domain input signal having a first spectrum at a 
first pre-determinable frequency or range of frequencies to 
25 a time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or range of 
frequencies, said system comprising 

means for identifying the first spectrum of the time or 
spatial domain input signal from which energy is to be 
30 transferred. 



4 
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means for specifying the second spectrum of the time or 
spatial domain output signal to which said energy is to be 
transferred, and 

means for calculating, using a frequency domain description 
5 of said output signal, for example, the output spectrum, 
expressed in terms of a frequency domain description of said 
input signal and coefficients of a time or spatial domain 
description of a generalised non-linear system, the 
coefficients of a time or spatial domain description of said 
10 generalised non-linear system in order to give effect to the 
energy transfer. 



A fourth aspect of the present invention provides a 
computer program product for designing a non-linear system 

15 for transferring energy from a time or spatial domain input 
signal having a first spectrum at a first pre-determinable 
frequency or range of frequencies to a time or spatial 
domain output signal having a second spectrum at a second 
pre-determinable frequency or range of frequencies, said 

20 computer program product comprising 

computer program code means for identifying the first 
spectrum of the time or spatial domain input signal from 
which energy is to be transferred, 

computer program code means for specifying the second 
25 spectrum of the time or spatial domain output signal to 
which said energy is to be transferred, and 

computer program code means for calculating, using a 
frequency domain description of said output signal, for 
example, the output spectrum, expressed in terms of a 
30 frequency domain description of said input signal and 
coefficients of a time or spatial domain description of a 
generalised non-linear system, the coefficients of a time or 
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spatial domain description of said generalised non-linear 
system in order to give effect to the energy transfer. 

A fifth aspect of the present invention provides a non- 
5 linear system which can transfer energy from a time or 
spatial domain input signal having a first spectrum at a 
first pre-determinable frequency or range of frequencies to 
a time or spatial domain output signal having a second 
spectrum at a second pre-detejrminable frequency or range of 
10 frequencies^ said system comprising 

means for identifying the first spectrum of the time or 
spatial domain input signal from which energy is to be 
transferred, 

means for specifying the second spectrum of the time or 
15 spatial domain output signal to which said energy is to be 
transferred, and 

means for giving effect to the energy transfer using 
coefficients of a time or spatial domain description of a 
generalised non-linear system, said coefficients having been 
20 calculated using a frequency domain description of said 
output signal, for example, the output spectrum, expressed 
in terms of a frequency domain description of said input 
signal and coefficients of a time or spatial domain 
description of a generalised non-linear system. 



Advantageously, the fifth embodiment allows processing 
for the determination of the coefficients to be performed 
off-line and merely incorporated into a non-linear system 
which uses the coefficients. 



25 



30 
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Embodiments of the present invention will be described 
by way of examples only, with reference to the accompanying 
in which: 

figure 1 shows the effect of traditional signal 
5 processing by, for example, a linear filters- 
figure 2 illustrates signal processing according to one 
aspect of the present inventions- 
figure 3 depicts a further example of the signal 
processing according to the present invention; 

10 figure 4 shows a non-linear system arranged to give 

effect to the energy transformation shown in figure 3; 

figure 5 illustrates a further energy transformation in 
which energy is distributed over a wider frequency band; 

figure 6 illustrates the power spectral densities for 
15 the input and output signals of a designed nonlinear system; 

figure 7 illustrates the power spectral densities for 
the input and output signals of another designed nonlinear 
system; 

figure 8 depicts schematically a non-linear system; 

20 figure 9 shows a digital implementation of the non- 

linear system shown in figure 8; 

figure 10 illustrates the power spectral densities for 
the input and output signals of the nonlinear system which 
was obtained using a design which improves the filtering 

25 effect shown in figure 6; 

figure 11 illustrates the power spectral densities for 
the input and output signals of the nonlinear system which 
was obtained using a design which improves the filtering 
effect shown in figure 7; 



7 
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figure 12 shows the time domain input and output of the 
nonlinear system with the frequency domain filtering effect 
shown in figure 10; 

figure 13 shows the time domain input and output of the 
nonlinear system with the frequency domain filtering effect 
shown in figure 11; 

figure 14 depicts the structure of a designed nonlinear 
system; 

figure 15 illustrates the frequency spectrum of a 
signal to be processed using the present invention; 

figure 16 illustrates the result of a Fast Fourier 
Transform of an input signal having the spectrum shown in 
figure 15; 

figure 17 illustrates the results of the n-dimensional 
(n= 2 and 3) convolution integration for the spectrum shown 
in figure 16; 

figure 18 illustrates the output magnitude frequency 
response of a designed nonlinear system to an input signal 
having the frequency spectrum shown in figure 15; 

figure 19 shows the frequency spectrum of a further 
signal to be processed using the present invention; 

figure 2 0 illustrates the output magnitude frequency 

response of the same nonlinear system as shown in figure 18 

to the further signal to be processed having the frequency 
spectrum shown in figure 19; 

figure 21 illustrates the frequency spectrum of a still 
further signal to be processed using the present invention; 

figure 22 illustrates the output magnitude frequency 
response of the same nonlinear system as in figure 18 to the 
still further signal to be processed having the frequency 
spectrum shown in figure 21; 

S 
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figure 23 illustrates the structure of another designed 

nonlinear systems- 
figure 2 4 shows the continuous time realisation of the 

discrete time system in figure 23; 

figure 25 shows a mechanical implementation of the 
continuous time system in figure 24; 

figure 26 illustrates the result of a Fast Fourier 
Transform of the further signal shown in figure 19; 

figure 27 illustrates the output magnitude frequency 
response of a further designed nonlinear system to the 
further signal shown in figure 19; 

figure 2 8 illustrates the output magnitude frequency 
response of a still further designed nonlinear system to the 
signal shown in figure 15; 

figure 2 9 illustrates the output magnitude frequency 
response of the nonlinear system shown in figure 28 to the 
further signal shown in figure 19; 

figure 30 illustrates a flow chart for designing a 
nonlinear system according to an embodiment; 

figure 31 illustrates a flow chart for designing a 
nonlinear system according to a further embodiment; 

figure 32 shows the structure of a nonlinear filter 
designed based on specifications for both the magnitude 
and phase of output frequency responses; 

figure 33 depicts the input and output magnitude 
frequency characteristics of a specific nonlinear filter 
designed based on specifications for both the magnitude 
and phase; 

9 
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figure 34 shows the phase angle of the spectrum 
^2(J^) in figure 32 in the specific design case shown in 
figure 33 which reflects the phase response 
characteristic determined by the design; 

5 

figure 35 shows the phase characteristics of the 
linear phase FIR filter in the specific design case shown 
in figure 33; 

10 figure 3 6 depicts the discrete time model of a 

nonlinear filter designed to focus energy from two 
different frequency bands into a single frequency band; 

figure 37 shows the spectrum of an input signal of 
15 the nonlinear filter in figure 36; 

figure 38 shows the frequency response of the 
nonlinear filter in figure 36 to the input in figure 37,, 
which indicates an energy focus effect of the nonlinear 
20 filter; 

figure 39 illustrates the block diagram of a spatial 
domain nonlinear filter; 



25 figure 40 depicts the power spectral densities of an 

input and the corresponding output of the spatial domain 
nonlinear filter in figure 39; 

figure 41 depicts an spatial domain input and 
30 corresponding output of the filter in figure 39; 

figure 42 shows a one dimensional image to be 
processed by the spatial domain nonlinear filter in 
figure 39; 



35 
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figure 4 3 shows the one dimensional image obtained 
by processing the image in figure 42 using the spatial 
domain nonlinear filter in figure 39; and 

5 figure 44 illustrates a data processing system upon 

which embodiments of the present invention can be 
implemented . 

Referring to figure 1, there is shown the principle of, 
10 for example, traditional low pass, high pass, and band pass 
filtering. Figure 1 shows the power spectrum of a signal 
100 both before and after processing. The energy of a 
signal 100 to be filtered comprises two parts, namely, a 
first part 102 for further processing or of interest and a 
15 second part 104 which is of no interest. Typically, the 
second part 104 of the signal is attenuated which results in 
a second signal 106, The second signal 106 comprises the 
original or a copy of the first part 102 and an attenuated 
portion or an attenuated version of the second part 108 . 

20 

Figure 2 illustrates the principle of signal processing 
according to one aspect of the present invention. Figure 2 
shows the power spectriom of a signal 200 both before and 
after processing. The signal comprises a first portion 202 

25 and a second portion 204. The first portion 202 of the 
signal 200 is of interest for further processing or output. 
Accordingly, as a result of the signal processing using the 
present invention, the first portion 202 is retained and the 
energy in the second portion 204 is translated to another 

30 frequency band 206. 
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: DETAILED DECRIPTION 



The theory and method underlying the present invention 
will now be described in general terms in steps (i) to (vi) . 

(i) Determine the frequency spectrum of a signal to be 
processed, including the range of frequencies of the signal, 

(i±) Specify the frequency spectrum of the output signal. 

(i±±) Deteinnaine the structure of a Nonlinear Auto-Regressive 
model with exogenous inputs (NARX model) to ensure that the 
energy transformation between different frequency bands and 
other design requirements, for example, specifications for 
magnitude and/or phase of the output spectrum over the 
required output frequency band can be met or realised. 

The general expression for a NARX model is given by 

N 

y(k) = Syn(k) (1) 

n=l 

where y^(k) is a 'NARX nth-order output' given by 

^ P p + q 

yn(k) = x s <=pq(ii'-"'ip.c,)riy(k-ii)n^(k-i,) 

p=01i,lp^q=l i = l i:=p+l 

(2) 

with 

p-l-q = n, 1^=1,. ..,K^, i = l,...,p + q, and 2^ = S'" S 

i,Ap*q-i ii=i ip^^^i 

is the maximxam lag and y ( .) / u(.), and c^^{ .) are the 
output, input, and model coefficients respectively. A 
specific instance of the NARX model such as 

y(k) = 0 .3u (k~l)+0 .ly(k-l)-0 .OZu (k-l) u (k-D'O .04u (k-2) u (k-1) 
-0. 06y (k-l)u (k-3)-0 .08y(k~2)y(k'-3) 

may be obtained from the general form (1) and (2) with 

12 



SUBSTITUTE SHEET (RULE 26) 



25 



30 



WO 99/45644 PCT/GB99/00550 

CQi (1)=0.3, CIO (1) =0.7, C02 (1/ l)=-0.02, cq2 ('2 , 1 ) =-0 . 0 4 , 
cii (l,3)=-0.06, C20 (2, 3)=-0.08, else Cpg(.)=0 

5 Simplified designs will be considered below where the 

NARX model with only input nonlinearities is employed. 
However, it will be appreciated by one skilled in the art 
that the present invention is not limited to use in relation 
to only input nonlinearities. The present invention can 

10 equally well be used in circiomstances of both non-linear 
outputs and non-linear inputs and outputs. Equally, the 
invention is not restricted to realisation as a NARX model. 
The invention may be realised using many alternative model 
forms either in discrete time or continuous time. Models 

15 such as the Hammers tein and Wiener model, or continuous time 
models, for example, a nonlinear differential equation model 
could be used or any other model, including discrete or 
continuous spatial models, that can be mapped into the 
frequency domain. However, for each of the models the main 

20 design principle is the same. 

The NARX model with only input nonlinearities is given 
by equation (1) where 



yn(k)= 



Jc^(l„-.,ljrju(k-l,) for n>2 

=1 i=l 
K K 

2;c,o{l,)y(k-l,) + 2;c,,{lJu(k-l,) for n=l 



K 



(31 



The structure of the NARX model (1) and (3) is defined 
by the values of K"^, n=l, . . . , N, and, for each n (an 

integer between 1 and N inclusive) , involves terms of the 
form 

n 

Co„(li/"-A„)n^('^-li) li=l,...,K„, i = l,...,n, 
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when n > 2 and 

Cio(li)y(k-lJ , Co,(l,)u(k-l,) , l,=l,...,Kj 
5 when n=l in the model. 

The parameter N in the non-linear model (1) and (3) is 
associated with the realisability of the model required to 
give effect to the energy transformation. The ability to be 
10 able to realise the energy transformation is determined from 
the relationship between the input and output frequencies or 
frequency ranges of non-linear systems. 



The structure parameters K^, n = l,,..,N, are associated 
15 with the extent to which specific design requirements such 
as the magnitude and/or phase of the output spectrum over 
the required output frequency band can be satisfied. These 
parameters are iteratively determined as part of the design. 
The model could initially be assumed to be of a simple form 
20 in terms of these parameters. However, if the initial 
choice of parameters does not produce a satisfactory design 
the parameters are progressively or gradually revised 
according to the energy transfer effect of the resulting 
non- linear system. 

25 

For systems described by the NARX model (1) and (3), 
the relationship between the input and output frequencies or 
frequency ranges is given by 



30 



35 



^Y-fY„U^Y„., (4) 

where f^ denotes the range of frequencies of the output, and 
f^^ and f^^_^ denote the ranges of frequencies produced by the 
Nth-order and (N-l) th-order nonlinearities , and 
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i -1 



k=0 
k=0 ■ 



when 
when 



nb 



(a + b) 

nb 
(a + b) 



na 



(a + b) 
na 



<1 



>1 



(a + b) 
n = N and 



N-1 



1 = 



+ 1 



where [.] relates to or means take the integer part, 
na 
(a + b) 

Ij, =[na-k(a+b) ,nb-k(a4-b)] for k = 0 - 1 , 
I^. = [0,nb-i*(a + b) ], 

and the frequencies of the signal to be processed are in the 
range defined by the interval [a,b]. 



10 Given [a,jb] and the required output frequency range fy, 

the smallest N for the NARX model (1) and (3) which can 
bring about the specified frequency domain energy 
transformation can be determined from equation 4. 



15 (iv) Map the NARX model with the structure given in (iii) 
into the frequency domain to yield the frequency domain 
description. The frequency domain description is given in 
terms of the Generalised Frequency Response Functions 
(GFRFs), H^( jw^,---, jw^) , n = l,...,N, which, after this mapping, 

20 are specified in terms of time or spatial domain model 
parameters . 



The mapping of the NARX model ( 1 ) and ( 3 ) between the 
time or spatial domain and the frequency domain is given by 

25 

1 ^ 

n = l,...,N (5) 
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Therefore the frequency domain properties of the system can 
be completely defined in terms of the parameters Cpg (.) of 
the time or spatial domain description of the system. 

5 (v) The output frequency response of the non-linear system 
(1) and (3) is given by 

N 

Y(jw) = 2;Y„(jw) (6) 



10 



where 



Yn(jw) = ^{^, jHn(jWi,...,jwjnu(jw,)da, (7) 

with 

J(.)do-„ 

15 denoting an integration over the nth-dimensional hyper-plane 
w^+,---,+w^ = w 

Based on this relationship, the parameters in 
(j^lf * ' ' rj^n^ ^^""^^ * ' - which, due to the mapping 

20 performed in (iv) , are the same parameters as those in the 
time or spatial domain model are deteirmined. This step 
enables the shape of the output frequency spectrum Y(jw) to 
be defined which in turn ensures that the spectrum 
approaches the specified output frequency spectrum as 
25 closely as possible. 

Different design specifications can lead to different 
implementations of corresponding designs. 



(v.l) Firstly, given knowledge of the input spectrum 

30 U(jw) and the required output spectrum (jw) . Substituting 
(5) and (7) into (6) yields 
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1 A l/ViT ^ 

il- 2-Cio(li) exp (-jwljS 

Ml (8) 



n 

J exp[-j (w,l,+,--.,+w„l„)]nU(jw.)dcr„ 



W, +,••--+ w„—w 



Equation (8) enables the parameters associated with the 
time or spatial domain NARX model, 

5 CQn di, ^ . . f In) f li=l f ' ' ' fKj^/ . ' ' ^1^=1^ , , . rK^r n=l ^ , . , 

and 

*^10 ^'^l) ' -^1^^' • • • / ^1, 

to be determined as follows to implement the required 
design: 



25 



1) Based upon the equation 

Y(jw) = Z T^t^ ^ Co.(ll^"vln)X 

J exp[- j (w,l,+,--v+w^l,)]nu( jw,)da„ 



w, +,---, + w_ =w 



i=l 



determine parameters 

^On (-^1' • • ^ f In) f ' ' • • / • • • >--^n~^' . . . , Kn^ n=l, . . . , 

using a least squares routine to make the right hand side of 
equation (9) approach the specified output spectrum as 
20 closely as possible. 

The first term on the right hand side of equation (8) 
1 



is omitted from (9). The omitted term represents linear 
output terms in the time or spatial domain realisation of 
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the nonlinear system and these may not be needed to achieve 
the design at this step, hence they are omitted from (9) . 

2) In order to augment the performance of a filter designed 
5 as above, it may be desirable to also design a suitable 
linear filter H{jw) to improve the approximation to Y* ( jpvt) 
obtained in 1) above such that 

n=l \^^} liAo^l 

J exp[-j {w,l,+,.--,+w„lj]nU(jw.)da-„ 

can achieve a better approximation to Y* (jw) , As part of 
10 this linear design the parameters ciqUi), ij. . =l,...,Ki^ 

which were omitted from (8) to get (9), can be obtained from 
the parameters of the linear filter. 

Design 1 hereafter illustrates the design of a non- 
15 linear system using this first case. 

(V.2) Secondly, given the input spectrum U(jw), and a 

specified bound for the magnitude of the required output 
spectrum, Y^* ( w) . 

20 A bound {w) for the magnitude of the output spectrum 

Y(jw^) of the NARX model (1) and (3) can be expressed as 

Y^w) = Z ,^^,(n-i) k( jw„---, jwj|' |u|*..-*|u ( jw)| (10) 

n 

25 according to the result in Billings, S.A. and Lang, Zi- 
Qiang, 1996, A bound for the magnitude characteristics of 
nonlinear output frequency response functions, Part 1: 
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Analysis and computation. Int. J. Control, Vol.65, pp309- 
328, where 

|u|*--'*|u ( jw)| 

n 

denotes the n-dimensional convolution integration for the 
magnitude characteristic of the input spectrum and 

K{jw,,---,jwj|^ 

represents a bound for the GFRFs magnitude 



with wi, . , . f Wjj satisfying the constraint pv?^i + . 
For the NARX model (1) and (3), 

K{j^i.--vjwj|^ 

can be evaluated as follows 



2 |c,,(l,,--,lj| (11) 



Y"(w) = 



1- S^iodi) ^^P (-jwlj 

Combining (10) and (11) yields 
1 



l-S^iodi) exp (-jwlj 



Z ,^ Li> c juh---Hu(jw)| 

n=l \ ' ' 



(12) 



where 

E|con(lH-.ln)| 
liA„-l 

Equation (11) enables 

Cn= Z|con(ll.-.ln)| 
liA„=l 

and cio (-^l)^ Ii=l , - . - , to be determined for shaping the 
bound Y^{w) for Y(jw) in order to make this bound approach 

Y^*(w) . The procedure which can be used to achieve this is 
as below. 
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1) Based upon the equation 

^'^^^) = Z ,^^,(n-i) Cju|*.-.*|u(jw)| (13) 

n 

use a least squares routine to make the right hand side of 
5 the above equation approach y'^^ {w) especially over the 
frequencies or frequency range to which the required energy 
transformation needs to be implemented. The coefficients C^^ 
n = i,---,N, in equation (13) must be constrained to be 
positive since C^^ is the result of the summation of the 
10 modulus of the coefficients cqjj (1 1 ^ . , , ^ Ijj) ^ 

1 J— 1 /-•••/ -fvn ^ ' • ' / -^n~"^ / • ' ' f -f^n . 

The first term on the right hand side of equation (12) 
1 

Ki 

15 is omitted from (13) . The omitted term represents linear 
output terms in the time or spatial domain realisation. 
These omitted terms may not be needed to achieve the design 
at this step, hence they are omitted from (13) . 

20 

2) If necessary, the approximation to Y^*(w) in 1) above 

can be supplemented using a linear filter with a magnitude 
25 characteristic |H(jw)| such that 

|H(jw)|i; Cjuh-*iu(jw)| 
n=i ( 2;r) . • 

n 

provides a better approximation to the specified bound and 
as a result cio(-^l)' li = l,.,.,Ki, which were omitted from 
30 (12) to get (13), can be obtained from the linear filter 
parameters . 
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Design 2 hereafter illustrates the detailed procedure 
of the above and several examples of the design. 

(V.3) Thirdly, there are many practical situations which 

5 should be dealt with on an individual basis. There follows 
examples of two such situations. 

(a) Referring to figure 3 there is shown an input 
signal having a spectrum U(jw) 300 which comprises a portion 
10 302 of interest, between frequencies a and e, and an portion 
304, between frequencies e and b, to be translated to 
another frequency range 306 defined by frequencies f and 2b 
while retaining the portion of interest 302 within the 
output signal 308. 

15 

In order to realise the energy transformation shown in 
figure 3, a non-linear system 400 can be constructed as 
illustrated in figure 4 which comprises means 402, 404 and 
40 6 for implementing the following components 

20 

Hi(jw), H2 ijwi, jw2) f and H(jw) . 

H^(jw) 402 and H{jw) 406 can be implemented readily using 
classical linear band pass filters, while ^2 ( j w^i , jiV2 ) 
can be constructed in the time or spatial domain as shown in 
25 the equation below 

y2(k)= X Co^d,, 12)0^(^-1.) 

with the parameters, Cosdi/lz) . li = l/"vK2,l2 = 1,---,K2 being 
determined to produce the signal yjCk) with a required 
30 frequency characteristic. Each of the components 402, 404 
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and 405 have corresponding frequency responses 408, 410 and 
412. Consequently, the whole non-linear system can be 
realised as a non-linear time or spatial domain filter as 

y(k) = H,(q-Mu(k) + H(q-') £ Co^C l,,!^)]^^ ( k - 

h.h=l i=l 

HiD^q ) ^d(^ ^ i=i 

where 

H.ijCq ) H^(q ) 

are the backward shift operator descriptions of the linear 
filters which have the frequency response functions Hi (jw) 

and H(jw), respectively. 

The equation for y (k) can be further written as 
H,i,(q'')H^(q'')y(k) = H,^(q-')H^(q-')u(k)+H„(q-MH,^{q-') £ ^02(1^,12)11^(^-1^ 

which is clearly a NARX model that can be described by 
equations (1) and (3) . 

Although the above shows energy transformation from a 
frequency band to a higher frequency band, energy can 
equally well be transferred from one frequency band to a 
lower frequency band. 

(b) If the objective of the energy transformation and 
hence the desired non-linear system is only to distribute 
energy of the signal to be processed over a wider frequency 
band without energy amplification then a simpler model may 
be sufficient. Referring to figure 5, there is shown 
schematically an input signal 500 in the frequency domain 
comprising energy between frequencies a and b. The desired 
output frequency spectrum 502 comprises two portions, a 
lower frequency portion 504 and a higher frequency portion 
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505. It will be appreciated that the energy of the input 
signal 500 is to be spread over the two frequency ranges 504 
and 50 6. 

5 A quadratic filter, for example, 

y(k)=a u2 (k) 

can be designed to redistribute energy of the original 
signal u (k) , with frequency components over the frequency 
10 band [a,b], to the new ranges [0,b-a] and [2a, 2b] without 
energy amplification provided an appropriate a is selected 
in the design process described above. 

(vi) The final step in the design process is to 

15 materially realise, that is, to physically realise the 
designed filter using appropriate software or hardware or 
combination of software and hardware. 

Although the present invention has been described above 
20 with reference to a NARX model, it will be appreciated that 
the present invention is not limited thereto. The method 
and realisation of non-linear systems having pre- 
determinable frequency or energy transfer characteristics 
can equally well be utilised using other forms of 
25 descriptions of non-linear systems or models where there 
exists a requirement to transform or transfer energy at one 
frequency or band of frequencies to energy at another 
frequency or band of frequencies . 

30 The present invention can be used to realise energy 

transformations over frequencies using a nonlinear system or 
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to modify energy transformation of an existing linear or 
non- linear system. 

Furthermore, the present invention can be applied or 
5 utilised in the field of design and realisation of 
electronic circuits or filters. The energy in the input 
signal at a particular frequency or band of frequencies is 
transferred to desired frequency bands. Similarly, in 
mechanical systems, the addition of non-linear mechanisms 
10 could transfer the energy of a vibration at an undesirable 
frequency to some other frequency. The present invention 
may also find application in the field of fluid mechanics, 
for example, in the effects of flow around objects {e.g. oil 
platform legs), noise in ducting and pipe flow systems. 

15 

Alternatively, the modifications which are required to 
be effected to a known linear or non-linear system, for 
example, a mechanical system, in order to bring about a 
particular frequency distribution of energy can be 
20 determined using the present invention and then the linear 
or non-linear system can be so modified. 

II . EXAMPLES 
25 II. 1 DESIGN 1 

Design No . 1 shows an example of how the energy of an 
input signal having predetermined frequency components can 
be transferred to other frequencies. 

30 

Consider digitally filtering an input signal 
u{t) = cost-hcos2t, (14) 
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wherein the sampling period is T = 1/100 s. 

The first step is to determine the frequency spectrum 
of the signal to be processed. The frequency spectrum of 
5 the signal to be processed comprises input frequencies 
w^^ = 1 and = 2 

The second step is to specify the output frequency 
characteristics of the filter. Two different filtering 
10 problems will be considered for this example. 



The specification for the first problem is to transfer 
the energy in u(t) to the output frequency w^^=0 and the 
specification for the second problem is to transfer the 
15 energy in u{t) to the output frequency w^^ = 4 . 

Since the sampling period is T = 1/100 s, the digital 
input to the filter is 
1 2 

u(k) = cos — k + cosj^k for k = 0,l,... , (15) 

20 the normalised input frequencies are therefore. 

w^i =1/100 and w^^ =2/100, 

and the required normalised output frequencies are 

w^^ = 0 for the first filtering problem and w^^ =4/100 for the 

second filtering problem, 

25 

For this example, the output frequencies produced by 
the nth-order nonlinearity of a non-linear filter are 
distributed uniformly in 
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fv = 



j = 0 



when 
when 



nb 



(a + b) 
nb 
( a + b) 



( a 
(a 



^ ]<1 
+ b) J 

+ b) J 



(16) 



1 2 ■ 

where ^ = ^0*0' ^^lOo"' ^ ^ ^ means to take the integer part 
of X, and 

na 



a. = 



(a + b) 



+ 1 



Ij =[na- j (a + b) , nb-j(a + b)] for j = 0, , 
I.. =[0,nb-r(a + b)], 

with the difference between any two neighbouring frequencies 

9 11 

being A = 



100 100 100 



10 



15 



20 



25 



that 



For n=l, it can be evaluated from the above equations 



lOO'lOO. 

and the corresponding output frequencies are 

_1 2_ 

lOO'lOO 

For n=2, it can be similarly obtained that 



(17) 



^ Lioo'loo. 



0,- 



100 100 



100 100 

and the corresponding output frequencies are 
12 3 4 



0, 



100 



0, 



lOO'lOO'lOO'lOO 



Therefore, a NARX model of up to second order 
nonlinearity will be sufficient to realise the energy 
transformations required by the filtering example thereby 
addressing the third step in the design process. 

According to the above analysis, select a NARX model of 
the form 



y (k)= 2]cio(ki)y (k-k:,) + X!Zco2(ki,k2)u(k-ki)u(k-k2) 



(18) 



ki=0k,=O 
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as the basic structure for the non-linear filter and take 
K2=l for simplicity. The parameters K^, c^q(1) c^q(K^) and 
Co2(k^,k2) f k^=0,l and k2=0,l, then need to be determined to 
completely specify the design. 

In order to derive the procedures for determining these 
parameters, consider the frequency domain characteristics of 
the filter model (18). According to J.C.Peyton Jones and 
S.A.Billings, Recursive Algorithm for Computing the 
Frequency Response of a Class of Non-linear Difference 
Equation Models, Int. J. Control, 1989, Vol.50, No. 5, 1925- 
1940, the generalised frequency response functions of this 
filter model are 



H,(jw) = 0 

1 1 

2 Z ^02( k„k2) exp[- j(w,k, + 



w 



(19) 




Hn( j^i/"*/jw„) = 0 for n>3 



The output frequency response of the filter (18) 
under the input 



1 , 2 , ^ A(w,) 

(k) = cos k + cos k= > ^^''^^ 

100 100 ,=_tl.o 2 



(20) 




can be 



written 



as 



N 




I EA(Wi^)A(WiJH2(jWi^,j^i^) for w->0 



Y(jw) = X!^n(jw) = Y2{jw) = 



n=l 
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where ±^,±2 e {- 2 ,-1 ,+1 ,+2} , Y(jw) is related to the output 
spectrum Y(jw) by 

r2Y( jw) for w >0 

[Y(jw) for w=0 

and, similarly, Yj^(jw) is related to the nth-order output 
spectrum Y^(jw) by 

r2Y,( jw) for w >0 

Y^(jw) = <^ 

I Y^(jw) for w = 0 

This is a. specific form of the general expression for output 
frequency responses of non-linear systems which are given by 

N 

Y(jw) = XY„(jw) (22) 

n-1 

where 

^n^^^^ = 7T^ K(jw,,...,jwjnU(jw,)dc7„ (23) 

with 

J{.)da-„ 

denoting an integration over the nth-dimensional hyper-plane 
Wi+,---,+w^ = w and N being the maximum order of the dominant 
system nonlinearities . 
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Substituting (19) into (21) yields 



1 1 



Y(jw) = i 2:A(w,^)A(w,^)^P^ 



1- Ecio(ki) exp[- jw kj'"=^'^=° 
k,=l 



11]^ r 
ZZco2(J<i'k2)- J^A^,^)A(vt^J exp[- j (w.^k, +w,^k2) 



1 1 



= H(jw)5]2^Co2(k„k2) f (w,k„k2) for w>0 



k,=0k2=0 



5 and 



f(w,k.,k,) 

Y(jw) = H(3w)2]Zco2(ki,k2) -i-^ for w=0 

k,=Okj=0 2 



where 



H( jw)= 



l-2c,o(k,)exp[-jw ki] 



(24) 



(25) 



(26) 



10 



f (w,k,,k2) = ^ XA(Wi^)A(w,^)exp[-j (Wi^k.+Wi^k^)] 



(27) 



15 



Moreover denoting 
i 1 



Y ( jw)=^ 



^^o^^ik^.k^) f {w,k,^,k^) for w>0 

for w = 0 



Jc,=0k2=0 
1 1 



gives 

Y(jw) = H(jw)Y(jw). w>0 



(28) 



(29) 



20 



In view of. the fact that H(jw) is the frequency 
response function of a classical linear filter and Y(jw) is 
a linear function of the filter parameters <^02^^if^2'> r 
k3^=0,l and kj^O,!, the procedures for determining the 
parameters of the non-linear filter with the given design 
requirements and the structure (18) are given as : 
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(a) From the design requirements, determine the desired 
output frequency characteristic Y*(jw) and choose the 
parameters CQ^i^if^z"^ f k^=0,l and k2=0,l, appropriately to 
make Y(jw) approximate y*(jw) as well as possible. 

5 

(b) Examine the filtering effect of 
1 1 

y(k)= 2] 2^02(^1^2) u (k-k2)u(k-k2) , (30) 

ki=0k2=0 

where CQ^ik^.k^) , k^=0,l , k2=0,l, are the results obtained 
in the procedure (a) . If the effect is acceptable, then 
10 choose H(jw)=l so that the filter parameters 

Cio( 1 ) = Cio(2 ) c,o(K,) = 0 ; 

otherwise design the classical linear filter H(jw) to make 
H(jw)y(jw) satisfy the requirements for the output frequency 
characteristics and obtain the filter parameters 
15 K^, Ciq(I) ,Cio(2 ) ,,..,Cj^q(K3^) at the same time. 

Therefore, in order to address the first problem to 
transfer energy from u(t), = 1 and \^^^=2, to the output 

frequency w^^ =0 , take 

20 Y*{jO) = l, Y*(jw,) = 0, Wi=Y^/ i-1.2,3,4. 

and for the second problem to transfer energy from u(t), 
w^3^ = 1 and w^2 = 2 , to the output frequency w^q=4, take 

25 The filter parameters c,^^(k^,k^) , k^=Q,l , k2=0,l, can then 
be determined through the group equations 
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k,=Ok,=0 



Re 



Y (j ) 

^100 



= ZSc,,(k„k,) f'*(--^,k„k,) 

k,=0)C2=0 XUU 



r— 1 v^v- f^(0,k,,k,) 



(31) 



k,=Ok,=0 



Im 



Y ( j ) 

-^100 



1 1 



= SZ=o2(ki,k,) f^(— ,k,,k,) 



by the least squares method. In (31) f'*( .) and f^( .) 
represent the real and imaginary parts of f ( . ) . 

Rewrite (31) as 



Y = X0 



(32) 
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Y = 



Re[Y*(jO)],---,Re 



Y (j ) 

100 



,Iin 



[Y*(jO)],. 



'flm 



Y (j ) 

100 



0 = [c^^{O, 0) ,Co2(l, 1) ,Co2(0, l) + Co2(l, 0)f 

f^O, 0, 0) f^O, 1, 1) f^O, 0, 1) 



X = 



f^'^TTT^T' 0,0), f^T-T-T/ 1,1), , 0,1) 



10 0 
f '(0, 0, 0) 



100 
f'(0, 1, 1) 



100 
f '(0, 0, 1) 



f'^TT^' 0,0), f '(r—r , 1,1), f'(T^' 0,1), 



100 



100 



100 



2 

2 c o s 



10 0 



2 c o s 



10 0 
2 

10 0 
J 

10 0 
4 

10 0 



2 1 

cos ^ + c o s 

10 0 10 0 

1 

cos 

10 0 
1 2 

cos + c o s 

10 0 10 0 

2 

cos 

10 0 
i 

2 sin 

S — 



10 0 



0 - 2 s i n 



0 - 2 s i n 



10 0 
2 

10 0 
_3 

10 0 

4 

10 0 



1 2 

sin sin 

10 0 10 0 

1 

-sin 

10 0 
1 2 

-sin sin 

10 0 10 0 

2 

— sin 

10 0 



(33) 



(34) 



(35) 



Notice that 0 can be written in the form of (34) due to the 
fact that 

f.(w, 0, 1)= f (w, 1, 0) . (36) 

Therefore, the filter parameters 0 for the first 
filtering problem to transfer energy from u(t), 
w^-L = 1 and w^2=2' "tihe output frequency w^^ = 0 can be 

obtained as 



0^ = (X'^Xj-^X^Yi 
where 



(37) 



1,0,---, 0,0,---, 0 



(38) 
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and the parameters 9 for the second filtering problem to 
transfer energy from u(t), ^ai = ^ ^a2 = 2 , to the output 

frequency w^q = 4 can be obtained as 



6i> = (X'*X)"'x''Y2 
where 



^2 = 



0,---,0 1 0,---,0 
The results of computing (37) and (39) are 





C'o2(0,0) 




1187 




C'o2(l, 1) 




1187 




_C'02(0, l) + C'o2(l, 0)_ 




-2373. 8_ 


and 










C'o2(0,0) 




'-I2O6.2' 


4 = 


C%2{1, 1) 




-1206.2 




C'o2(0, 1) + C'02(1, 0) 




2413.2 



(39) 



(40) 



(41) 



(42) 



respectively . 



Under the input of (20) , the power spectral densities 
of the input and output of the non-linear filter 
1 1 



yi(k) = X!Z^''>2(k,,k3)u(k-ki)u(k-k2) 

ki=0k2=0 

= c^2(0,0)u'(k) + c'o2(ia)u'(k - 1) + [c'o2(l,0) + c'o2(0,l)]u(k - l)u(k) 



(43) 



which is initially designed to address the first filtering 
problem above, are shown in figure 6, and the power spectral 
20 densities of the input and output of the non-linear filter 
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1 1 



= c'o2(0, 0)u'(k) + c'o2(l, l)u'(k-l) + [c'o2(l / 0) + c'o2(0, l)]u(k-l)u(k) 

(44) 

which is initially designed to address the second filtering 
5 problem above, are shown in figure 7. 

Further improvement to the performances of the filters 
(43) and (44) may be possible. Therefore, a linear filter, 
H(jw), is designed in order to further improve the filter 
10 performances . 



To improve the performance of the filter (43) for the 
first filtering problem, H(jw) is designed to be a fifth 
order low-pass type 1 Chebyshev filter with cut off 
frequency 0.5 rad/sec and 0.5 dB of ripple in the pass-band 
to ensure a satisfactory frequency response over the pass 
band. The result is 

b ^(l) + b^(2) z-^4---f-b^(6) 
l-a'(2) z"' a'(6) z"' 

^.,2 0. 0174 + 0. 08 53z"' + 0 . 1812z-^ + 0 . 1652z"^ + 0 . 0933z"' + 0 . 0159z"^ 



H,(jw)=- 



= 10" 



l + 4.9941z"^ -9.9765z'^ +9.9648z"^ -4.97 66z"'' +0.9942z"^ 

(45) 

To improve the performance of the filter (44) for the 
second problem, H(jw) is designed to be a fifth order high- 
pass type 1 Chebyshev filter with cut off frequency 3.9 
rad/sec and 0.5 dB of ripple in the pass band for the same 
25 purpose as in the first problem case. The result is 

b'(l) + b='(2) z-'+--+b'(6) 



H, ( j w ) = • , 
^ l-a'(2)z-' a'(6)--^ 



z 



0 . 9218 


-4. 608 8z'^+9.2175z"^ 


-9.2175z"^ + 4.6088z"'' 


-0.9218Z"' 


1 + 4 


. 8381z"^ -9. 3635z"^ +9 


. 0613z"^ -4 .3846z"' +0 


. 848 6z~^ 



(46) 



30 The purpose of the additional linear filter is to 

attenuate unwanted frequency components in the outputs of 
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the filters (43) and (44) for the above two filtering 
problems respectively to make the output of the additional 
filter satisfy the corresponding design requirement. 

It will be appreciated that the filter parameters 
K^, c^od) ,---,Cio(K^) , associated with the expression 




in H{jw) can be obtained as a result 



of dividing the denominator of H(jw) by the numerator of 

H(jw). The specific H(jw) is given by H,(jw) andH^Cjw) in 

10 equations (45) and (46) for the two filtering problems 
respectively . 



15 



20 



The general description for the nonlinear filters 
designed as above is 



y(k)-£c,o(k,)y(k-kJ = 



l-Z'=io(ki)q'''' y(k) = H(q-^)y(k) 



"2^ "2^ ^2 ^2 

==ZZCo2(ki,k2)u(k-k,)u(k-k2)= 5]Z^02(^^l/k2) ^T'^T'^W (47) 

where 

_ K, \ 

H(q"') = l-Zc,o(k,)q-"^ , (48) 

and q"S qi^ and q'^ denote the backward shift operators. 
Another expression in the time domain for (47) is 

y(k) = H-^(q-M[ 22co2(k,,k2)q-''>q-'== ]u(k) (49) 

kj=0k2=0 



25 Therefore, embodiments of the non-linear filters which have 
been designed can be realised in a manner as shown in figure 
8 or more specifically as shown figure 9. 

It will be appreciated that in figure 8 the two 
30 components 800 and 802 are represented by 

G 1 ( q 1 S q^ M = i] i c ^ J k, , k,)q;'^ q"^^ 

ki=0k2=0 
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and 

1-a (2)q-' a (6)q-' 

For the first filtering problem in this example 

with Co2(ki,k2) 's given by (41) and 

g ^ b^(l) + b^(2)q-^+--fb^(6)q-^ 

' l-a'(2)q-' a'(6)q-' 

with b^(i) 's and a^i) 's given by (45) . 



For the second filtering problem in this example 
10 G,(q-\q-^)= X;Sc^„,(k„k2)q,-'''q-''' 

with 602(^1, kj) given by (42) and 

^ , , _ ( 1 ) + ( 2 ) q-^ +. • -^b^ ( 6 ) q-^ 
^ fq ) . 

l-a'(2)q-' a^6)q-' 

with b^(i) 's and a^(i) 's given by (46) . 

15 Referring to figures 10 and 11, there is shown the 

filtering effects of the filters designed to address the two 
filtering problems. Figure 10 shows the power spectrum 
densities of the output and input of the nonlinear filter 
which is finally obtained for the first filtering problem. 

20 Figure 11 shows the power spectrum densities of the output 
and input of the nonlinear filter which is finally obtained 
for the second filtering problem. The effects of the 
filtering in the time domain of the two filters are shown in 
figures 12 and 13. All the responses indicate that the 

25 filters substantially satisfy the design requirements. 
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II. 2 DESIGN 2 

Design No . 2 illustrates the detailed procedure and 
several examples for designing a non-linear system for 
5 transferring the energy of a signal from a first 
predeterminable frequency or range of frequencies to a 
second predeterminable frequency or range of frequencies 
such that the output frequency response of the nonlinear 
system so designed is within specified bounds. 



II . 2 . 1 Detailed Procedure 

(1) Given u(t), the signal to be processed in the time or 
spatial domain, the frequency band [c,d] to which the 
energy of u(t) is to be transferred, and the user 
specified bound Y®*{w) for the output spectrum Y{jw) over 
[c,d] for the design. 

(2) Sample the time or spatial domain signal u(t) with 
sampling interval T to yield a discrete series {u(k)}and 
perform a Fast Fourier Transform (FFT) on the series to 
compute the spectrum U(jw) of u(t) as 



and M is the length of the data used to perform the FFT. 
M is taken as an even number for convenience. 

(3) Evaluate the range [a,b] of frequencies in u{t) as 



10 




where U^[j (•)] is the result of the FFT operating on {u(k)} 



2n 

b = 1 

MT 
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M 



and 1^ is an integer such that 

5 (4) The relationship between the bound of the output 
spectrum Y^(w) , the coefficients of the NARX model 

y(k) = 2]yn(k) 



10 



y„(k) = 



z-z c,„(i„-,i„)n^(k-ii) for ^^^2 

il=l i^t=l i=l 
Xcio(li)y {Jc-l,) + X^oi(li)u(k-l,) for n^l 



and the spectrum U(jw) is given by 



1 - S c,^( 1,) exp (-jwi^; 

1 



1,=1 1„ = 1 



Z .0^!'^n-l) Cjuh...^|U( jw] 



|U|*..-*|U ( jw)| 



where 



15 



= Z*'-ZKn(ll.->ln)l' n = N,,...,N 

are parameters associated with the NARX model parameters 

Kn/Co^(li.-'-/ln) ' li =lr-vK„. "-/In K„, for n = No,.'',N, 

|u|*...*|u ( j w)| denotes the n-dimensional convolution 

n 

integration for the magnitude |u(jw)|of the spectrum 
U(jw), which is defined by 
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|uh--*|u( jw) |= £-f^|u ( jw,)|...|u [ j (w-w,- .-.-^-w^.,) ]|dw,...dw„., 

n 

and = 1 when the NARX model involves nonlinear terms 
from order 1 to N . 



10 



15 



20 



25 



Based upon this expression, the structure 
parameters N and Nq and the NARX model parameters are 
determined as below. 



(i) Evaluate 



i=l 



Li=o 



when 
when 



nb 



(a-hb) 
nb 



na 



(a + b) 

na 
(a + b) J 



<1 
> 1 



1 = 



+ 1 



(a + b) 

where [ x ] denotes the integer part of x, 

na 
(a + b) 

I- =[na-i (a + b) ,nb-i (a + b)], for i = 0,--vi*-l 

I.. =[0,nb-i*(a + b)] 
for n=l,2, . . . until a value of n is reached such that 
part of the specified output frequency range [c,d] 



falls into f„ 



of Nq. 



This value of n is used as the value 



(ii) Evaluate 

for n=2,3, . . . until a value of n is reached such 
that the frequency range [c,d] falls completely 
within the corresponding fy . This value of n is 
taken as the value of N . 

(iii) Calculate 

lu|*...*|u(jw)l 
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to yield 

|u|*...*|u(j27r i/MT)|, for i = 0,--vM/2 

n 

using the algorithm 

|uh- (j 2n i/MT)| = Tg[i + f l]n ^ = 0 

n 

|u (0 ) ,"-,U [n (M - 1 )| = Conv|[u ( 0) ,"-,U (M - 1 )],--•, [3 (0 ) ,---,U (M - 1 )] |, 

' 27t( M V 

j i +1 

M V 2 J 



i = 0 , 1 ,"-,M - 1 



for n = No,---,N, where Conv ( . ) denotes the 

convolution operation and U(.) and U ( •) represent 
the intermediate results of this algorithm. 

(iv) Based on the (i^-i^+1) equations 



where i^ = round 



cMT 
2n 



i^ = round 



dMTl 
. 2;r J 



and round ( x ) 



means to take the integer nearest to use a least 
squares routine to compute 

, n = No,.-vN, 

under the constraint that the results must be 
positive, and then select the NARX model parameters 

for n = NQ,--vN 

under the constraints on the summation of the 
modulus of the coefficients given by 

Cn = S"-i;|con(li.-.lJ, n = N,,-.,N 

ll=l 1„=1 

(v) If necessary, design a classical linear filter, for 
example, a band pass filter which ideally allows the 
frequency response to be unity over the frequency 
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band [c,d] and zero beyond to yield the linear 
frequency characteristic 

1 

and therefore determine the parameters associated 
with K^, Cjod^), li=l,--vKj. Otherwise, all of the 
parameters of c^q( .) can be taken as zero to yield a 
model having no regression terms associated with the 
output . 

(5) Construct a NARX model as shown in figure (14) using the 
results obtained in the above (iv) and (v) . The 
nonlinear system illustrated in figure 14 comprises a 
nonlinear part 1400 and a linear part 1402. It will be 
appreciated that in figure 14 

X c^(lxr-vl.)flu(k-i,)=i;-..|; cjl,,-vl,)nu(k-l,) 

II. 2. 2 Three Specific Examples 
Example 1 

This example illustrates a further implementation or 
design of a nonlinear system following the above detailed 
procedure . 
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(1) The signal to be processed is given by 
, , ^ 1 sin a t -sin 6 t 

with a = 3,3, /3=^1, and = 1 . 6 . The frequency spectrum 
U(jw) of the signal is shown in figure 15 indicating 
5 that the real input frequency range is [l,3.3]. The 

requirement for the design is to transfer energy of the 
original signal to the frequency band [c,d] = [ 5 . 6 , 7 , 6] 
with the bound on the output spectrum magnitude 
specified to be Y^*(w) = 1.6 over this frequency band. 

10 (2) Sample u(t) with sampling interval T = 0.01 sec to produce 

1 sinakT-sin y5 kT ^ ^ ^ 1 sin 3 . 3 x 0. Olk-sin 0 . 01 k 

u(k) = 2iy[^- = 2x 1, 6x — X — 

2;r kT 27r 0.01k 

k = -1999,"-,0,---,2000 
and perform a Fast Fourier Transform (FFT) on this 
15 series (M=4000) to compute 

U.(j^l)=U,(j^l) 1 = -1999,...,0,...,2000 

and then to yield 

the result of which, in the nonnegative frequency range, 
20 is shown in figure 16. 

Notice the difference between the real spectrum of 
u(t) in figure 15 and the computed spectrum in figure 
16. The differences are due to the errors caused by the 
25 FFT operation. The design should and will be performed • 

based upon the computed spectrum to lead to more 
practical results . 
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10 



(3) Evaluate the frequency range [a,b] of u(t) from the 
computed spectrum giving 

a = 0.6283 b = 3,7699 

(4) Design the system structure and parameters 

(i) Determination of . 

Clearly, the output frequency range contributed by 
the linear part when the input frequencies are within 

[a,b] = [0. 6283, 3.7699] is 

f^^ =[a,b] = [0.6283, 3.7699] 

The frequency range f^^ produced by the second 

order nonlinearity in this case .is obtained as 
follows . 

Since n = 2. 



20 



nb 



15 (a+b) 



na 



2x3.7699 



(a+b) J (3.7699+0.6283) 



2x0.6283 



(3,7699+0.6283) 



and 



So 





na 




i* = 




+ 1 = 




Ja+b). 





=1.7143-[ 0.2857 ] = 1.7143-0>1 

+ 1 = [ 0.2857 ] + l = 0+l = l 



2x0.6283 



3.7699+0.6283 



^ =LX =] li^-i (a+tj)] Ll[0,ri^i*(a+b)] 

=[na,r±^yo,rb-(a+b)]=[2x0.6283, 2x3.769^|Jo, 2x3.7699-(0.6283+3.7699)] 

= [l.2566, 7.5398]U[0, 3.1416] = [0, 7.5398] 

thereby obtained contains part of the specified 

output frequency range [c,d] = [5.6, 7.6]. So, A''^ is 
determined to be Ng = 2 . 

(ii) Determination of N . 
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Evaluating fy = [Jf^ ^ f or n=2 yields 

f.|n=2=f.,U^. =[0' 7,5398]U[0.6283, 3.7699] = [o, 7.5398] 



To evaluate fy|n=3 = tf^U^z ' calculate f^^ first. I 
this case 



n 



nb 



na 



3x3-7699 



3x0.6283 



(a+b) L (a+b) J (3.7699+0.6283) [ (3.7699+0.6283) 

=2.5714-[ 0.4286 ]=2.5714-0>1 

3x0.6283 





na 


+ 1 = 


i* = 






.(a+b). 





3.7699+0.6283 



+ 1 = [ 0.4286 ] + l = 0 + l = l 



So 



=Uli =IoUli =[na,nb][J[0,nb-(a+b)] 

i=0 

= [3x0.6283, 3x3.7699]|J[0, 3x3. 7699-(0. 6283 + 3 .7699)] 
= [1.8848, 11.3097]U[0, 6.9115] = [0, 11.3097] 
Therefore 

f,|„^3 =f,^Ufy. =[0' 7.5398]1J[0, 11.3097] = [0, 11.3097] 

^¥^=3 thereby obtained includes the whole specified 

output frequency range [5.6, 7.6]. N is therefore 
determined to be N = 3 . 

(iii) Calculate |uh-**|u ( j w)j for n = = 2 and n = N = 3 , 

n 

respectively to yield 

M^-^U(j2;ri/MljHl^^-iu(j27ri/4000x0.0lj i=0,- - >400(/2, n=2 ard 3 
^ ^ ' V ^ , 

n n 

The results are shovm in figure 17, 

(iv) Based on the (i^-i^+1) equations 

Y^*(7r i/2000x0.01) = 1.6 = -^C2 |u|*|u(j ;r i/2000 x 0 . 0 1)| 
+ -^^C3 |u|*|u|*|u(j;r i/2000 xO.Ol)| 



with 
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10 



15 



= r ound[cMT/2;r] = r ound[5 . 6x4000x0, 01/2x;t] = 36 
= round[d]yiT/2;r] = round[7 . 6 x 4000 x 0 . 01/2 x ;r] = 48 
use a least squares routine under the constraint of 
nonnegative solutions to compute C2 and C3 , that is, 
to determine C, and C3 , under the constraints of 



C2 > 0 and C3 > 0 , 
expression 



to 



minimises 



the 



following 



^ r I 1 

^1.6-— q |i:j^i^j;r 1/2000x0. 0l)|-^C3 |c|Hu(*|u(j^ 1/2000x0. Ol)| 



The results obtained are 



C2 =0 



and 



C3 =3. 8367 



(v) Design, optionally, a linear Butterworth band pass 
filter to attenuate the frequency components beyond 
the frequency band [c, d] = [5.6, 7.6] to yield the 
linear frequency characteristic 



10-'(0 . 098 6-0 . 1972q-' + 0 • 098 6q-') 
1-3 . 9 633q"^ + 5 . 8988q"^ -3. 9076q~^ + 0 . 9721q" 



q = e^ 
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(5) Construct a NARX model as shown in figure 14 with 

No=2, N = 3, K2-K3=l, Co2(l.l) = 0, Co3( 1 , 1 , 1 ) = 3 . 8 3 67 , 
that is^ 

Z Z---Z^(1i'"A)E^(K-1J =Z£<^(vi2)]l^(k-lJ+ 

n=^, 1^1=1 a„=l i=i li=ll2=l i=l 

K, K, 3 

SEE ^3(ii/i..i.)ni(3c-ij 

11 2 

=zz co2(ii.i2)n^(k-ii)+ 

li=ll2=l i=l 

111 3 

zzz co3(ii.i2.i3)n^(k-i,) 

Il«ll2=ll3«i i=l 

=qj2(l,l)u^k-l)+qj3(l,l,l)u'(k-l)=0xu'(k-l)+3.8367u'(k-l)=3.^^^ 

and 

lCr'(0 . 0986-0 . 1972q-' +0 , 0986q^) 
1-3 . 9633q-' +5 . 8988q"^ -3 . 9076q"' +0 . 9721q^ 
which determines the parameters associated with the NARX 
10 model parameters and c^gdj , 1^=1, •-•^K^. 

The output frequency response under the given input is 
shown in figure 18 indicating that the energy has been 
transferred to the specified frequency band [c, d] = [5.6, 7.6] 
15 with the magnitude of the response below the specified bound 
1.6. 

The frequency response of the above design is examined 
for other input signals below. 



1,=.! 
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In the first case, consider 



u(t) = 



2M 



2cos 



(a + b) t 
2 



-cosbt-cosa t 



;r{b-a) t^ 



where a, b, and are defined as above. The frequency 

spectrum of the signal is shown in figure 19, The frequency 



5 range is clearly the same as that of the signal having the 
spectrum given in figure 15 and the magnitude of the 
spectrum also satisfies the condition that 



This implies that the frequency response of the designed 
10 system to this u(t) should in theory also transfer energy 



magnitude frequency response being less than Y^*(w)=:1.6 over 



response and indicates that the actual result is consistent 
15 with the theoretical predictions. 

For the second case, u(t) was taken as a random 
process with the frequency spectrum given in figure 21. The 
frequency spectrum is substantially within the frequency 

20 range [l, 3,3] with a magnitude of less than 1.6. Therefore, 
the same conclusion should apply for the output magnitude 
frequency response of the designed system to this random 
input. Figure 22 shows this response and indicates that the 
energy is transferred to a new frequency band of 

25 substantially [5.6, 7.6]. Note that the magnitude of the 
output spectrum over this frequency band is well below the 
specified bound Y®*{w) = l,6. This is because of the effect of 
the attenuation which is due to intrakernel interference of 
the nonlinear mechanism. 

30 The nonlinear filter which has been designed above can 

be represented by the block diagram shown in figure 23. This 



|u(jw)|<M^=1.6 



into the frequency band [c, d] = [5.6, 7.6] 



with the output 



this frequency band. 



Figure 2 0 shows this frequency 
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could be realised electronically by following the approach 
used to realise Design 1 in figure 9. However, for some 
applications the design may need to be realised in 
continuous time. Using the bilinear transformation 

_ 1 + (t/2)s 
1-(t/2)s 

where T=0.01 is the sampling interval, q"^ is the delay 
operator, and s is the Laplace Transform operator and 
substituting in the discrete time design in figure 23 
provides the equivalent continuous time system shown in 
figure 24. Simulating the system in figure 24 produces 
almost identical responses as the discrete time system 
response of figure 18. 

The system in figure 24 could for example be realised 
mechanically as illustrated schematically in figure 25 where 
the cubic device is either a material which exhibits a cubic 
response or is implemented as an actuator which takes u^^Ct) 
as input and produces an actuation output UjCt) which is 
proportional to the cubic power of the input. 

One possible application of this design would be in 
vibration isolation. For example it may be required to 
transfer energy from an input frequency range [1,3.3] to the 
frequency range [5.6,7.6], The new design could be used to 
achieve this effect. 

Other much more complicated designs for vibration 
isolation can be achieved based on the present invention. 
The design procedure would be exactly as described above. but 
the realisation would involve the synthesis of dampers, 
damping materials or actuators with the nonlinear dynamic 
characteristics specified by the designs. 

Example 2 
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Example 2 shows an application of the design above to 
the attenuation of signal energy over unwanted frequency 
bands using designed nonlinear effects. 

When using linear structures for the attenuation of 
5 signal energy in a physical system, the attenuated energy is 
usually absorbed by devices such as dampers in mechanical 
systems and resistors in electronic circuits and transformed 
to other energy forms such as thermal energy. This may lead 
to undesirable effects and measures, such as using radiating 

10 devices, sometimes have to be taken to compensate for these 
effects. When a nonlinear system is employed, instead of 
attenuating the signal energy directly as in the linear 
case, the signal energy at a frequency of interest can be 
spread over a wider frequency band and attenuated by means 

15 of the counteractions between different terms which compose 
the output spectriam. This means that, to a certain extent, a 
nonlinear design for signal energy attenuation can reduce 
the requirements for using energy absorption devices and is 
of great benefit in practical applications. 

20 

Another important application would be, for example, in 
the design of the foundations or the modification of the 
characteristics of buildings and structures which are in 
earthquake zones. The objective in such an application would 

25 be to design materials or actuators, tuned as required for 
each structure, which transfer the damaging input energy 
from an earthquake to another more acceptable frequency band 
or spread the energy to be within an acceptable bound over a 
desired frequency range. Spreading the energy using the 

30 present design should produce significant reductions in 
earthquake damage . 



(1) Design a nonlinear system to attenuate the energy of the 
signal 
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2M. 



;T(a-b) 



2cos 



fa-fb ^ 



cosa t - cosh t 



where a = 3.3, b = l, M^ = 1.6. 
(2) The spectrum of the signal is evaluated. This is 
implemented by sampling the signal with sampling 
interval T = 0,01 sec to produce 



u{k) = 



2M. 



n{a- p) ( kT f _ 



2cosl — - — kT -cosor kT-cos>9 kT 



(3) 



k = -1999,---, 0,-", 2000 
and perform a Fast Fourier Transform (FFT) on this 
series. The result of the FFT is shown in figure 26. 

Evaluate the range [a,b] of frequencies in the signal 
based on the computed spectrum. The evaluation gives 



a = l,0996 b=:3.1416 

This is because the computed spectrum indicates 

|u(jw)|>0,05 for w €[l. 0996, 3,1416] 

and 

|u(jw)|<0.05 for other w 

(4) Assume that, as part of the design, output frequency 

range is [c,d] = [o, 7.3] and the required bound over this 
frequency band is Y^*(w) = l, 

(5) Following the same steps for the design of the system 
structure and parameters as in example 1 above, 
Nq, N, and C^, n = NQ,---,N were determined as follows. 



(i) No 



is determined to be Nq = 1 . This is because 



f,^ =[a,b]-[l. 0996, 3.1416]g[0, 7.3] = [c, d] 
where f^^ denotes the output frequency band contributed 

by the system linear part, and part of the selected 
output frequency range falls into the linear output 
frequency range f^^ . 
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15 



20 



(ii) N 

In this case, it can be obtained using a=1.0990 and 
b=3.1416 that 

f^^ =[2.1992, 6.2832]|J[0, 2.0420] 
f,^ =[0, 9.4248] 

Thus, if the maximum nonlinear order is taken to be 2, 
then the output frequency range of the system is 

f.|n=2=f^Uf.. =[0' 6.2832] 
If the maximum order is taken to be 3, then the output 

frequency range 

f.ln=3=f.3Uf.. =[0' 9.4248] 

Clearly, fr|n=3 contains the whole selected output 

frequency range [c,d] = [O, 7 . 3] , N is therefore determined 
to be N=3. 

(iii) C^, n = l,2,3. 

C^, n = l,2,3. are determined to minimize the 
following expression 

gi-cJu{j2;zi/Mrj-^q |u|*|u{j2;riyiytr)|--^C3 luHuj*|u(j2;r i/MT)! 



under the constraints of C^>0, i = l,2,3. 

Notice that for this specific example, 
M = 4000 



i^ = round 



ij = round 



CMT 

dMT l 
. 2^ J 



j = rounc 



= round 



0 X 4000 X 0 . 01 



7 . 6 X 4000 X Q . 01 
2fc 



= 0 



= 48 



The solution to this minimisation problem is 
C,=0, C2=2.1932, C3 = 6.0550 
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(6) Select the NARX model parameters 

^Qn^^if'-'flJ , 1^ =l,-'-,K^r",l^ -1,'-,K^, for 

n = No,---/N , 

based on the results obtained in (5) as 

5 K, = = K3 = 2 

cbi(.) = 0, 

q,2(l,l) = 1.1932, qj2(1^2) = c^{2,l) = 0, c^^{2,2)=--l, 

Cij3(l,l,l) = 3.0550, qj3(2,2,2) = -3, and other 0^3(0 's are zero 
and consequently construct a NARX model as 

y(^)=E S"*Z^on(ii."vijn^(k-i,)=S^oi(ii)u(k-i,)+ 

n=rNo ii=l la=l i=l 

K2 2 K K K 3 

co2(i„i^)nu(k-i,)+ 1;!;!; co3(i„i„i3)riu.(k-i,) 

ll=U2=l i=l lj=il2=H3=l 

11=02=1 i=l ii=lX,=li,=l i=X 

= q32(i'l)^'(k-l)+q,2(2,2)u'(k-2)+q,3(l,l,l)u'(k-l)+q33(2,2,2)u'(k-2) 
= 1. 1932u^(k-l)-u^(k-2) + 3 . 0550u^(k-l)-3u^(k-2) 

Notice that the selected NARX model parameters 
15 satisfy the relationship 

E|coi(lx)| = 0 = C, 

ii=i 

Z 2^2(11.12)1=^2(1/ 1)1+^2(2,2)1 =|l -19321+1-1] =2.1932=C2 

l£SWVl..J3Hq«B(l.l,lHq3(2,2,2j H3.055CH-:^ =6.0550=0, 

20 and the different signs selected for 

Co^d.l) and Co2(2,2) 

and for 
Co3(l,l,l) and Co3(2,2,2) 

are such as to give effect to the intra-kernel and 
25 inter-kernel interferences, which is to attenuate the 

energy of the input signal. 
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The frequency domain response of the constructed model 
to the sampled series of the input signal is shown in figure 
27 . It can be seen that excellent energy attenuation has 
been achieved by the designed system. It will be 

5 appreciated that the input energy in figure 2 6 has been 
spread over the designed frequency band by the nonlinear 
filter . 



Example 3 



10 



This example shows another application of the present 
invention to the attenuation of signal energy over unwanted 
frequency bands using designed nonlinear effects. The 
example also illustrates the effect of the same design on a 
15 different input signal to demonstrate the effectiveness in 
energy attenuation of the designed system in different 
circumstances . 

(1) Design a nonlinear system to attenuate the energy of the 
signal 

/V ^.,1 sin a t -sin (5 t 

20 u t =2M„ 

" 2;r t 

with a = 3. 3, ^ = 1, andM^=1.6. 

(2) The spectrum of the signal is evaluated by sampling the 
signal with sampling interval T=0.01 sec to produce 

' , ^ 1 sin cd^-sin kT ^ 1 sin 3.3x0. Olk-sin 0.01 k 

u(k) = 2M^- ~ =2x1. 6x — X 

2;r kr 2jt 0.01k 

25 k = -1999,---, 0,---, 2000 

and then performing a Fast Fourier Transform (FFT) on 
the obtained time series. The result of the FFT is the 
same as that shown in figure 15. 
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(3) Evaluating the range [a,b] of frequencies in the signal 
based on the computed spectrum gives 

a = 0. 6283 b = 3.7699 

as the computed spectriom indicates 

|u(jw)|>0.05 for w.e[0,6283, 3.7699] 

and 

|u(jw)|<0.05 for other w. 

(4) Assume that, as part of the design, it is desired that 
the output frequency range is [c,d] = [O, 10 . 3] and the 

required bound over this frequency range is Y^*(w)=l. 

(5) Following the same steps for the design of the system 
structure and parameters as in example 1 above, the 
parameters Nq, N, and C^, n = NQ,-",N are determined as 
follows . 

(i) No 

Nq is determined to be Nq = 1 due to 

=[a,b] = [0.6283, 3.7699]e[0, 10.3] = [c, d] , 
which indicates that part of the selected output 
frequency range falls into the linear output 
frequency range f^^ . 

(ii) N 

In this example, N is obtained using a=0.6283 and 
b=3 .7699 . Therefore, 
f,^-[0, 7.5398] 

fy^ =[0, 11.3097] . 

Hence, if the maximum nonlinear order is taken to be 
2, the output frequency range of the system is 

f.|n.2=fv,U^: =[0' 7,5398]. 

If the maximum order is taken to be 3, the output 
frequency range 

f.|n=3=f..Uf.. =[0' 11-3097]. 
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15 



20 



contains the whole selected output frequency 

range [c,d] = [o , 7 . 6] , . N, for the design, is determined 
to be N=3. 

(iii) C^, n = l,2, 3 

C^, n = l,2,3 are determined in such a manner as 
to minimise 

2l-q|u(j2;ri/Mr(-^C, \x^^d,j2m/m)\--^C, |[^Hu(-|l^j2;r i/Mr)| 

under the constraints of C.>0, i = l,2,3. In this case 
also M=4000 but 



i^ = round 



i^ = round 



cMT 

dMT 
2^ 



fO X 4000 X 0 . 01 

= round 

L 27r 



= round 



= 0 

3 X 4000 X 0 . 01" 



= 66 



The solution to this minimisation problem is 
C,=0, 0^=0. 2928, C3 =0.9763 
(6) The NARX model parameters 

^nf<=^oJ^xr",lJr 1, =l,---,K„,-"/l„ =1/---,K„, for n = No,-.-,N, 
are selected based on the results obtained in (5) as 

q,:(.)=0 

q^(l,l) =0.1928, c^ar2)=-0.1, q,2(2,l) = 0, qj2(2,2) = 0, 
q,3(l, 1,1) =0.6763, q33(l,2,2) = -0.3, and other qj3(.) 's are zero, 
and consequently a NARX model is constructed 

y(k)= Z cjl„-,ljfl^(k-lj 

n=No 1^=1 i=i 

• =Co2(l, l)u^(k-l) + Co2(l,2)u(k-l)u(k-2) + Co3(l,l,l)u^{k-l) 
+ Co3(l,2,2)u(k-l)u'(k-2) 

=0.1928u^(k-l)-0.1u(k-l)u(k-2)+0.e763u^(k-l)-0.3u(k-l)u^(k-2) 
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The selected NARX model parameters satisfy the 
relationships 

i;ki(ij|=o=q 

ii=i 

Z2^2(li'i2)| =1^)2(1' 1)1+19)2(1' 2^ =|0 •1926! +1-0 -l] =0.2928=C, 

ii=i i,=a 

and the different signs selected for 

Co2(l/ 1) and Co2(l / 2) 
and for 

Co3( 1,1,1) and Co3(l,2,2) 
10 are in order to give effect to the intra-kernel and 

inter-kernel interferences, which is to attenuate the 

energy of the input signal. 

The frequency response of the constructed model to the 
15 input signal specified above is shown in figure 28 which 
indicates that the required energy attenuation has been 
realised. 

The frequency response of the design above to the input 
20 signal in example 2 gives the result shown in figure 29. 
The nonlinear system design above clearly also works for the 
signal in example 2 in energy attenuation although the model 
was not specially designed for this signal. This is 

reasonable since the magnitude of the spectrum of the signal 
25 in example 2 is less than the magnitude of the spectrum of 
the signal in this example over almost all of the input 
frequency band and over the other frequency bands the 
magnitudes of the spectra of the two signals are all zero. 
This illustrates that the above design is effective not only 
30 for the input based on which the design is implemented but 
also for other inputs with magnitude frequency 
characteristics less than the magnitude of the spectrum of 
the considered input. 
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III . FLOWCHARTS 

Referring to figure 30, there is shown a flow chart 
3000 depicting the steps of an embodiment of the present 
5 invention in which the output spectrum is specified over a 
given range of output frequencies. 

The signal to be processed and the desired frequency 
response of the non-linear system to be designed are input 
10 via steps 3002 to 3006. 

At step 3002, a digital input signal {u(k)} and its 
sampling interval T are to be provided. The range of output 
frequencies [c,d] over which the energy of the input signal 
15 is to be transformed is given in step 3004. The output 
frequency range is specified using beginning and end 
frequencies c and d respectively. The distribution of the 
energy over the output frequency range [c,d] is requested or 
specified in step 3006. 

20 

The frequency characteristics of the input signal are 
determined at steps 3008 and 3010. More particularly, the 
frequency components of the digitised input signal, {u (k) } , 
are calculated using a Fast Fourier Transform at step 3008. 
25 The range of frequency components contained within the input 
signal is determined from the FFT at step 3010. 

Referring to steps 3012 to 3024, the orders of the 

nonlinearities required to realise a desired nonlinear 

30 system and, hence, the energy transformation, are 
calculated. 
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Variables n and Nq are set to one and zero respectively 
at step 3012. The output frequency components resulting from 
the nth-order of system nonlinearities are determined at 
steps 3014 and 3016. 

5 

A determination is made at step 3018 as to whether or 
not Nq=0, that is to say, whether or not the smallest order 
of system nonlinearities which makes a contribution to the 
desired energy transformation has been determined. If Nq is 

10 equal to zero, a deteintiination is made at step 3020 as to 
whether or not part of the specified output frequency range 
falls within fyn- "^^^ determination is negative, 

processing continues at step 3026 at which the value of n is 
increased by one. However, if the determination is 

15 positive, the value of is set to equal n at step 3022 and 
processing continues at step 3026. 

If Nq is not equal to zero, then a determination is 
made at step 302 4 as to whether or not the specified output 
20 frequency range lies completely within fy = fyn "^-^Yn-l* 

the determination is positive, processing continues with 
step 3027. However, if the determination is negative, the 
value of n is incremented by one at step 3026 and processing 
continues at step 3014. 

25 

Steps 3028 to 3034 determine the values of the lags of 
the nonlinear model for all values of n = Nq, No+i,...,N and 

determine the parameters of the nonlinear system to be 
designed. 

30 Step 3036 represents a step for fine tuning of the 

output or frequency response of the designed non-linear 
system. 
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Steps 3038 and 3040 determine whether or not the 
resulting frequency response is sufficient to meet the 
design requirements. If the frequency response of the 
designed filter is acceptable, the discrete version of the 
5 filter is output at step 3042. However, if the designed 
non-linear system does not have an acceptable frequency 
response, then step 3044 increments the value of by one 

for all n= Nq, No+i,..-,N and steps 3030 to 3040 are 
iteratively repeated . 

10 

The determination in step 3040 as to whether or not the 
designed filter is acceptable is made, for example, by 
calculating the difference between the required output 
spectrum and the real output spectrum of the designed system 

15 over the frequency band [c,d]. If the modulus of the 
difference is below a predeterminable threshold over all the 
output frequency band [c,d], then the designed nonlinear 
system may be deemed to be acceptable. However, . if the 
modulus of the difference is greater than the threshold at 

20 any frequency over [c,d], then the design is refined. 

Referring to figure 31, there is shown a flow chart 
3100 for implementing computer code according to a second 
embodiment of the present invention. 

25 

The signal to be processed and the specification for 
the desired bound to be imposed upon the output signal 
frequency characteristics is input at steps 3102 to 3105. 

30 At step 3102, a digitised input signal and its sampling 

interval T are to be provided. The range of output 
frequencies [c,d] over which the energy of the input signal 
is to be transformed is specified in step 3104, A bound 
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Y^* (w) for the distribution of the energy over the output 
frequency band [c,d] is input at step 3106. 

The frequency characteristics of the input signal are 
determined at step 3108 and 3110. More particularly, the 
frequency components of the digitised input signal, { u (k) } , 
are calculated using a Fast Fourier Transform at step 3108. 
The range of frequency components contained within the input 
signal is determined from the FFT at step 3110. 

Referring to steps 3112 to 3124, the orders of the 
nonlinearities required to realise a desired nonlinear 
system and, hence, the energy transformation, are 
calculated. 

Variables n and Nq are set to one and zero respectively 
at step 3112. The frequency components resulting from the 
nth-order of system nonlinearities are determined at steps 
3114 and 3116. 

A determination is made at step 3118 as to whether or 
not Nq = 0, that is to say, whether or not the smallest 
order of system nonlinearities which contributes to the 
desired energy transformation has been determined. If Nq is 
equal to zero, a determination is made at step 3120 as to 
whether or not part of the specified output frequency range 
falls within f Yn • determination is negative, 

processing continues at step 3126 in which the value of n is 
increased by one. However, if the determination is 

positive, the value of Nq is set to equal n at step 3122 and 
processing continues at step 3126. 

60 



SUBSTITUTE SHEET (RULE 26) 



wo 99/45644 




PCT/GB99/00550 



If, at step 3118, it is determined that Nq is not equal 
to zero, then a determination is made at step 3124 as to 
whether or not the specified output frequency range lies 
completely within fy fyn ^ ^Yn-l- ^^e determination is 

5 positive, processing continues with step 3127, However, if 
the determination is negative, the value of n is incremented 
by one at step 3126 and processing continues at step 3114. 

Steps 3128 to 3134 determine the values Cj^ for all 

10 values of n = Nq/ N0+I/-../N, which represents the summation 

of the modulus of the values of the parameters of the 
nonlinear system. 

At step 3136 the parameters of a NARX model are 
selected given the constraints imposed upon the non-linear 

15 system obtained in steps 3128 to 3134. A determination is 
made via steps 3138 and 3140 as to whether or not the 
frequency response beyond the frequency band [c,d] of the 
designed filter is acceptable. If the frequency response is 
acceptable, the filter design parameters are output at step 

20 3142. However, if the filter characteristics are not 
acceptable at frequencies outside the frequency range [c,d], 
then a conventional filter, H(q"^), is designed, at step 
314 4, in order to reduce or obviate the frequency response 
of the designed filter outside the range of frequencies 

25 [c,d]. Finally, the design is completed at step 3146 by 
combining the designed non-linear and linear filters, if 
any . 



The deteirmination in step 3140 as to whether or not the 
30 frequency response of the designed filter is acceptable 
outside the range of frequencies [c,d] is made, for example, 
by comparing the modulus of the frequency response beyond 
the frequency range [c,d] with a predeterminable threshold. 
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If the modulus is below the threshold over ail the frequency 
range outside [c,d], then the designed nonlinear system may 
be deemed to be acceptable. However, if the modulus is 
greater than the threshold at any frequency beyond [c,d], 
the design is refined. 
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IV. THREE MORE DESIGNS AND EXAMPLES 

IV. 1 DESIGN OF NONLINEAR FILTERS WITH SPECIFICATIONS FOR 
BOTH THE MAGNITUDE AND PHASE OF OUTPUT FREQUENCY RESPONSES 

5 

It will be appreciated that the basic principles of the 
present invention can be applied to the design of nonlinear 
filters based on specifications for both the magnitude and 
phase of output frequency responses. The ability to 

10 modulate phase as well as or instead of magnitude is of 
particular importance in telecommunication applications. 

Consider a filtering problem such that given an input 
spectrum U(jw) over a frequency band [a,b] and a desired 

output spectrum Y\jw) over another frequency band [c,d], a 
15 filter is required to be designed so that the output 
frequency response Yij^) can match the desired spectrum 

as closely as possible in teirms of both magnitude and 
phase characteristics. The basic principles in Section I can 
be directly applied to realise the design of such a filter. 
20 The procedure to be followed is described below. 

First, a nonlinear filter 

y,(t) = N[uit)] 

is designed using the basic principles described in Steps 
(i) — (iv) and (v.l) Part 1 in Section I to produce a 
25 frequency response Yiijw) such that Yiijw) can match Y*(jw) 
over the specified output frequency range [c,d] as closely 
as possible in terms. of both the magnitude and phase. 

Secondly, if required, design a linear filter with a 
frequency response function H,(jw) such that, ideally, 

30 H,(jw)^^^^j^ wG[c,d] 

Y,(j^) 
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This can be used to improve Y^{jw) , the result obtained from 
the nonlinear design, so that the frequency response of 
the corresponding output 

Y2Uw) = H,(jw)Y,(jw) 

5 provides a better match to the desired spectrum Y^^ijw) . 

Thirdly, if required, design a linear phase FIR (Finite 
Impulse Response) band pass filter H2(jw) with the ideal 
magnitude frequency characteristic 



\h,Uw)\ = 




w G[c,d] 
otherwise 



10 and a linear phase over the frequency range. Then construct 
the designed filter using linear filters H^(jw) and H^Cjpx?-) 
and nonlinear filter N[u(t)] as shown in figure 32 to yield 
the output frequency response 

15 The second and third steps above follow the design 

principle described in Step (v.l) Part 2 in Section I so as 

to augment the performance of the nonlinear filter 
Yiit) ~ N[u(t)] designed in the first step. 

Ideally, Y(jw) = H2{jw)Y^_(jw) implies that 

20 \y(M^\ ==\H,Um^U^i-\y\jw^ w^[cd\ 

[0 otherwise 

and 

ZY(jw) - ZH^{jw)'\-ZY2U^) = k^w-^ZY\jw) 
due to the linear phase characteristic of H^Cjw) where is 
a coefficient which is a function of the order of H^ij^) . 
25 This indicates that the output frequency response of the 
designed filter is ideally the same as the desired response 
over the output frequency range [c,d] except for a linear 
phase difference between the real and the desired phase 
response. This is in fact an unavoidable phenomenon if band 
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Figure 35 shows the phase angle of the applied 
linear phase filter H^{jw) . 

It can be observed from the above that a nonlinear 
filter designed using the basic principles of the present 
5 invention can produce an output frequency response which 
satisfies a design specification in terms of both the 
magnitude and phase . 

The discrete time model description of the designed 
filter is given below where the nonlinear part of the 
10 model is the discrete time model description of the 
nonlinear filter 



and the linear part of the model is the discrete time 
model description of the linear filter, the frequency 
15 response function of which is 



Nonlinear part of the model: 
20 y(k)= +(3.123e+06) u(k-l) + (-2 . 244e+07 ) u(k-2) +(6.602e+07) 



y,(t) = N[u(t)] 



H(,jw) = H^(,jw)H,SD^) 



u(k-3) 



+ (-1 . 027e+08) u(k-4) +(8.961e+07) u(k-5) +(-4.18e+07) 



u (k-6) 



+(8.176e+06) u(k-7) +(4.239e+09) u(k-l)u(k-l) 



25 



+(-1.772e+10) u(k-l)u(k-2) +(5.897e+09) u(k-l)u(k-3) 
+(7.213e+09) u(k-l)u(k-4) +(-3.66e+09) u(k-l)u(k-5) 
+(-1.869e+08) u(k-l)u(k-6) +(6.736e+07) u(k-l)u(k-7) 
+(-6.777e+08) u(k-2)u(k-2) +(8.003e+10) u(k-2)u(k-3) 
+(-7.57e+10) u(k-2)u(k-4) +(6.023e+09) u(k-2)u(k-5) 
+(8.648e+09) u(k-2)u(k-6) + (-4 . 3 65e+08 ) u(k-2)u(k-7) 
+ (-8 . 365e+10) u(k-3)u(k-3) +(2.268e+10) u(k-3)u(k-4) 
+(1.005e+ll) u(k-3)u{k-5) + (-3 . 795e+10) u(k-3)u(k-6) 
+(-2.516e+09) u(k-3)u(k-7) +(1.054e+ll) u(k-4)u(k-4) 
+(-1.97e+ll) u(k-4)u(k-5) +(l.l09e+10) u(k-4)u(k-6) 



30 
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+(1.907e+10) u(k-4)u(k-7) +(2.352e+10) u(k-5)u{k-5) 
+(8.173e+10) u(k-5)u(k-6) + ( -3 . 318e+10 ) u(k-5)u(k-7) 
+(-4.218e+10) u(k-6)u(k-6) +(2.042e+10) u(k-6)u(k-7) 
+(-1.658e+09) u(k-7)u(k-7) ^ 

Linear part of the model: 



10 



15 



20 



25 



30 



35 



y(k)=b(l)u(k)+b(2)u(k-l)+, . . . , +b (m) u ( k-m) -a { 1 ) y ( k-1 ) . . 
a (n) y (k-n) 

with n=2 and m=303 where 

[a (1) , , a (n) ] = 

-1 .8 62158716053 98 0 . 95715116734794 
and 

[b(l) , , b (m) ] = 

l.Oe-03 * 
Columns 1 through 4 

0.00000003328202 0.00000021432147 0.00000080368128 
0.00000227704991 

Columns 5 through 8 

0.00000537879856 



0.00001113499106 



0.00008548478679 



0. 00024351205626 



0.00003571014847 

Columns 9 through 12 

0.00005704570668 
0 . 00016141668710 

Columns 13 through 16 

0.00020406372559 
0.00028246765302 

Columns 17 through 20 

0.0002 6339508245 0 . 0002 057274 4049 
-0. 00005450964923 
Columns 21 through 2 4 

-0. 0002 618125042 9 -0 . 000513818 00480 
-0. 00108568381434 
Columns 25 through 2 8 

-0. 00135214 637123 -0 . 00155851835832 
-0. 00162976908378 
Columns 2 9 through 32 



0 . 00002079938253 



0 . 00012079738459 



0 . 00027258483896 



0 . 00010137844659 



-0.00079603373996 



-0.00166438952001 
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-0.00141994747987 -0.00101090816448 -0.00039475153522 

0.00041552374840 

Columns 33 through 36 

0.00138259540242 0.00244344402372 0.00351062196838 

5 0.00447655873267 

Columns 37 through 4 0 

0.00522101718800 0.00552147075427 0.00556579258464 

0 .00496627276545 

Coliamns 41 through 4 4 
10 .0.00377365278712 0.00198 963381599 -0.00032378155160 
-0.00303968275897 
Columns 45 through 48 

-0.00596769040539 -0.008 8 6230677 97 6 -0.01143843394688 
-0 . 01339358387193 
15 Columns 49 through 52 

-0.01443551408603 -0.01431325012233 -0.01284878916639 
-0 . 00996641286493 
Columns 53 through 56 

-0.00571602453045 -0.00028726506096 0.00598841256370 
20 0.01264965179911 

Columns 57 through 60 

0.01913182739200 0.02480771068003 0.02903991461588 

0 . 03124135106508 

Columns 61 through 64 

25 0.03093862246516 0.02783233132424 0.0218 47 8 68 63383 

0 .01317044094614 

Colijmns 65 through 68 

0.00225896524855 -0.01016501602446 -0.02315515892661 
-0. 03560414378179 
30 Columns 69 through 72 

-0.04633007891956 -0.05417943987784 -0.05813822823951 
-0 . 05744137376473 
Columns 73 through 7 6 

-0.05166958483783 -0.04082303832663 -0.02536256454818 
35 -0.00621129546254 

Columns 77 through 8 0 

0.01528703636150 0.03745307958913 0.05839584942467 

0 . 07616729599331 
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10 



15 



20 



25 



30 



35 



Columns 81 through 84 

0 . 08 8 935377 01557 0 , 0951607 84 83677 

0.08424478625381 



0 .09376093809767 



Columns 8 5 through 88 

0.06680350720636 0.04234527006789 0.01246669428018 

-0.02064087852648 

Columns 8 9 through 92 

-0.05434613862312 -0.085783374 61088 -0.11209316760456 

-0.13067820746590 

Columns 93 through 96 

-0.13945154456062 -0.13705390892471 -0.12301823047738 

-0.09786316009190 

Columns 97 through 100 

-0.06310299056015 -0.02116847625143 



0. 07098889554710 
Columns 101 through 104 



0. 11362138287669 
0.18477606135825 



0 . 14887942949652 



Coliamns 105 through 108 

0. 1813498 4027062 0 . 1628 4 982 658 307 

0. 08570350882143 



0 . 02475895107438 



0.17344622849594 



0 - 13023456167487 



Columns 109 through 112 

0.03255064587105 -0.02508361142904 -0.08254318260252 
-0 . 13504269472806 
Columns 113 through 116 

-0.17807650738416 -0.20781448143133 -0.22144885867665 

-0.21746018536730 

Columns 117 through 120 

-0.19577680884658 -0.15781155123986 -0.10636991902139 
-0 . 04543567185132 
Colximns 121 through 124 

0. 02014 93057 6277 0.08508 4147 98056 

0. 19217566873157 

Columns 125 through 128 

0. 2254244124 9983 0. 24097 17 1607 669 



0.14405304899329 



0.23743645368538 



0.21500740460840 

Columns 129 through 132 

0. 1754325 6598 964 0 . 121874 9924 77 61 



0 . 05864768640035 
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-0.00914923229050 
Columns 133 through 136 

-0.07605725937466 -0.13671663683833 -0.18631209098366 
-0.22096869885059 
5 Columns 137 through 140 

-0.23806415895179 -0.23 643147511610 -0.21643603250516 

-0.17992231311682 

Columns 141 through 144 

-0.13003701658436 -0.07094606546323 -0.00747191949549 

10 0.05531598034943 

Columns 145 through 14 8 

0.11252149812387 0.15981666793597 0.19377819943088 

0.21213874239435 

Columns 14 9 through 152 

15 0.21393429471004 0.19953878954962 0.17058700949905 

0. 12979660938955 

Columns 153 through 15 6 

0.08070832589097 0.02736960335274 -0.02 6008 82 98 9955 
-0.07539394440027 
20 Columns 157 through 160 

-0.11724060322409 -0.14874679880287 -0.16803574024719 
-0. 17425382855726 
Columns 161 through 164 

-0.16758132308460 -0.14916007746448 -0.12094944070741 
25 -0.08552673296083 

Columns 155 through 168 

-0.04585223365109 -0.00502017035384 0.03398322 634818 

0. 06849574758055 

Columns 169 through 172 

30 0.09635411698109 0.11601725927068 0.12663256434628 

0. 12804494388799 

Columns 173 through 17 6 

0.12075316710299 0.10582180562013 0.08475985344594 

0. 05937858562038 

35 Columns 177 through 18 0 

0.03164147930605 0.00351816237707 -0.02314740959980 

-0.04674674074459 

Columns 181 through 184 
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-0.06598945884170 -0.07995651694641 -0.08812411270222 
-0. 09036059631364 
Columns 185 through 188 

-0.08690017680683 -0.07829817708583 -0.06537295841595 
5 -0.04913952047510 

Columns 18 9 through 192 

-0.03073930999045 -0.01137008308799 0.00778109613199 

0.02559935810198 

Columns 193 through 196 

10 0.04109626424759 0.05345544314259 0.06206856688305 

0.06656076607150 

Columns 197 through 200 

0.06680468888031 0.06292254225138 0.05527568217195 

0.04444169619049 

15 Columns 201 through 204 

0.03117946545587 0.01638338990942 0.00102875697948 

-0. 01388895099522 

Columns 205 through 208 

-0.02741724189660 -0.03870846951445 -0.04707559699096 
20 -0.05203784736309 

Columns 20 9 through 212 

-0.05335260858275 -0.05103084289941 -0.04533443158680 

-0 . 03675527219878 

Columns 213 through 216 

25 -0.02597741626758 -0.01382495500742 -0.00119959263130 

0. 01098723154804 

Coliamns 217 through 220 

0.02188229583806 0.03075319656401 0.03704059636285 

0. 04039599306472 

30 Columns 221 through 224 

0.04070212260137 0.03807460516092 0.03284505620907 

0. 02552746332356 

Columns 225 through 228 

0.01677103005022 0.00730378144712 -0.00212808430644 
35 -0.01081988478115 

Columns 229 through 232 

-0.01816221252030 -0.02368447649585 -0.02708498340864 
-0 . 02824561802129 
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Coliimns 233 through 236 

-0.02723060984541 -0.02427029083261 -0.01973202360355 

-0.01408149907660 

Columns 237 through 24 0 

5 -0.00783828203366 -0.00152977531128 0.00435233046871 
0.00938952420923 

Columns 241 through 244 

0. 01326503243370 0.01578162483958 0.0168 6791315533 

0.01657347478885 

10 Columns 245 through 24 8 

0.01505411203293 0.01254934451451 0.00935476560966 

0. 00579214028759 

Columns 249 through 252 

0.00218008037598 -0.00119216975076 -0.00408581800286 
15 -0.00632711824010 

Columns 253 through 256 

-0.00781483560799 -0.00852066157236 -0.00848329376098 
-0.00779740923517 
Coliomns 257 through 260 
20 -0.00659902170090 -0.00504908820003 -0.00331691201488 
-0.00156482692685 
Columns 2 61 through 2 64 

0.00006463440339 0.00145824730714 0.00253783620427 

0.00326220436512 

25 Columns 2 65 through 2 68 

0.00362556118310 0.00365304002511 0.0033 940 9451670 

0.00291465798950 

Columns 269 through 272 

0.00228895177653 0.00159174429109 0.0008 9171267597 

30 0.00024636506918 

Columns 273 through 27 6 

-0.00030123128510 -0.00072387 920912 -0.00100998481168 
-0. 00116188515783 
Columns 277 through 280 
35 -0.00119315515203 -0.00112532046764 -0.00098439594641 
-0.00079761965498 
Columns 2 81 through 284 

-0.00059067279083 -0.00038 55777567 6 -0.00019936402666 
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10 



15 



WO 99/45644 

-0 . 00004349568836 

Columns 285 through 288 

0.00007602486265 0 
0. 00022267635728 

Columns 289 through 292 



00015801946885 



0. 00021739726997 
0. 00013324716819 



0.00019645030203 



Columns 293 through 296 



0.00010088214188 
0. 00003090906954 



0.00007223559481 



Columns 297 through 300 



0.00001822614592 
0. 00000207399807 



0 . 00000988083116 



Columns 301 through 303 

0. 0000007 42 0 9218 0. 0000002 0073804 
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0. 00020513691086 



0.00016652267489 



0. 00004878213424 



0 . 00000483478363 



0. 00000003164528 



TABLE 1 
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IV. 2 DESIGN OF NONLINEAR FILTERS TO FOCUS ENERGY FROM 
DIFFERENT FREQUENCY BANDS INTO A SINGLE FREQUENCY BAND 

The above embodiments of the present invention 
5 illustrate how the principles underlying the invention can 
be used to transfer energy from one frequency or frequency 
range to another frequency or frequency range • However^ it 
will be appreciated that the present invention can equally 
well be utilised to focus energy from a given frequency 
10 band or given frequency bands into a single, and 
preferably, narrower frequency band. 

The problem to be addressed in this design is that 
given an input spectrum which possesses nonzero magnitude 
characteristics over two different frequency bands [a,,jb,] 

15 and [a2yb2], where a3 > jbj , a filter is required to be designed 
to focus energy from the two different input frequency 
bands into a single output frequency band [c,d], where 
c>jb, ,d<a2, with the spectrum of the filter output over the 
output frequency range [c,d] satisfying certain 

20 specifications. The specifications to be satisfied can be 
specified in teimis of magnitude, phase or both. 

The only difference of this design compared to the 
designs in Section I is that the determination of the 
maximum order N of the filter nonlinearities should not 

25 only ensure that the specified output frequency band [c,d] 
is covered by the filter's output frequency range, but also 
should ensure that the output energy over [c,d] is derived 
from the input frequency range [a,,jbj and the input 
frequency range [a^.b^] ^ Therefore, the principles for this 

30 design are exactly the same as those described in Section I 
except for step (iii) which should follow the new principle 
above . 
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Consider an example where a nonlinear filter is to be 
designed to focus the energy from two different frequency 
bands [3,4] and [10,11] of an input signal 



271 t In 



where = 4, = 3, jb^ = 1 1, a. = 10, = 0.03 , into a single frequency 

band around frequency f=7 with the single output frequency 

band as small as possible and the output magnitude 
characteristic at frequency f=7 being 3M^ = 0.09 . 

It can be shown that when subject to the excitation of 
an input with the frequency spectrum over two different 
frequency bands [a,,jbi] and {^a^.h^^, the output frequency 
ranges of a nonlinear system contributed by the system' s 
2nd order nonlinearity is 

\la^:ih,\ [2a, ^,], [a, +a3,i), +it,], [0,i?, -aju[0,jb, -a^] and \e^~h,,h.^-a,\ 

15 Substituting jb, = 4, a, = 3, jb^ = 1 1, a^ = 10 into a-, - iDj and - 
yields 

[a, -h,,h,_ -ai] = [6,8] 
Clearly this output frequency range covers frequency f=7 
and the system output energy over this frequency band is 
20 derived from the input frequency range [^i = [3,4] and the 

input frequency range [a2,b-;] = [lO^l] . Therefore, the maximum 
order of system nonlinearities can be taken to be N=2 for 
this specific design. 

Moreover, following the same design principles as 
25 described in Section I, a discrete time nonlinear filter 
under the sampling period r^ = l/50s can be obtained for this 
specific design. Figure 36 illustrates schematically such 
a filter. 
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Figures 37 and 3 8 show the input and output frequency 
spectra of this filter respectively and clearly indicate 
that the input energy over the two different input 
frequency bands [3,4] and [10,11] have been, as required, 
5 successfully focused into a very small band around 
frequency f =7 . 

IV. 3 DESIGN OF SPATIAL DOMA.IN NONLINEAR FILTERS 

Although the above embodiments have been described 
10 with reference to the design of non-linear time-domain 
systems, it will be appreciated that the present invention 
is equally applicable to non-linear spatial-domain systems • 

In the one dimensional case, the present invention can 
be directly applied and the only difference of the one- 
15 dimensional spatial domain case over the above time-domain 
embodiments is that the time-domain variables t and k in 
the continuous and discrete time filter equations become 
the continuous and discrete spatial variables X and in 
the spatial domain filter equations. 

20 Consider a one dimensional spatial domain filtering 

problem that is the same as the second problem addressed in 
Section II, 1 except that the time variables t and k become 
the spatial variables X and in the present design. 

Figure 39 shows the block diagram of the spatial 
25 domain nonlinear filter designed in this case. It can be 
seen that the spatial-domain filter is exactly the same as 
that shown in figure 9 which is the block diagram of the 
corresponding time domain nonlinear filter designed in 
Section II, 1 except that the time variables t, k, T, and 
30 s (second) in figure 9 have been changed to the spatial 
variables X, X^,AX, and m (meter) in figure 39. 

Figures 4 0 and 41 show the effects of spatial domain 
nonlinear filter. Again, except for the time variables 
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being changed to spatial variables, the figures are the 
same as figures 11 and 13 which are the results of the 
corresponding time domain nonlinear filter. 

The filtering effect in the time domain of the second 
5 nonlinear filter designed in Section II. 1 indicates an 
energy transfer from the input frequencies and w^^. = 2 

to the output frequency w-^^ = 4 . Figure 40 indicates the same 
effect but in the spatial domain. In this case, the spatial 
domain nonlinear filter transfers energy from the input 
10 spatial frequencies w^, = 1 and w^2=^ to the output spatial 
frequency w^q = 4 . 

When this one dimensional spatial domain nonlinear 
filter is applied for one dimensional image processing, 
where figure 41 represents the intensities over the spatial 

15 variable X of the images before and after the processing, 
the energy transformation effect of the spatial domain 
nonlinear filter can be further illustrated by figures 42 
and 43. Figure 42 shows the one dimensional input image 
with energy located in spatial frequencies w^^=l and ^^2~^ 

20 and figure 43 shows the one dimensional output image with 
energy located in another spatial frequency w^o = 4 due to 
the effect of nonlinear filtering. 



Although the present invention has been described 
above with reference to a one dimensional case in the 
spatial domain, it will be appreciated that the present 
invention is not limited thereto. The present invention can 
equally well be applied to the design and implementation of 
non-linear systems for realising transfer of energy from 
first m-dimensional spatial domain frequencies to specified 
second n-dimensional spatial domain frequencies, where m 
and n are greater than one. 
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It will be appreciated that an application of the 
transfer of energy from first m-dimensional spatial domain 
frequencies to specified second n-dimensional spatial 
domain frequencies would be digital image processing or 
5 filtering. In such a case m=n=2 . The present invention 
can be designed to produce filters that operate upon 
digital images. The filters can be designed to perform 
numerous different functions such as, for example, image 
compression or filtering to remove noise or vary the colour 
10 space of an image. 

The phrase "transfer of energy" includes, without 
limitation, the processing of a first signal, specified in 
the time domain or spatial domain, to produce second signal 
having predeterminable characteristics including a 
15 specified energy distribution. 

Furthermore, it will be appreciated that the terms 
"frequency range" and "range of frequencies" includes a 
group of frequencies or frequency group, A group of 
frequencies comprises a plurality of frequencies spatially 
20 distributed within an n— dimensional space or over a 
subspace . 

It will be appreciated by those skilled in the art 
that the phrase "specified spectrum'" relates to the 
specification of at least one of either the magnitude and 
25 phase of a signal or frequency components of a signal and 
is equally applicable to situations in which only the 
magnitude is specified and to situations in which both the 
magnitude and phase are specified. 

It will be appreciated that the input and output 
30 signals of any or all of the various embodiments of the 
present invention may be time or spatial domain continuous 
or discrete input or output signals. 
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Referring to figure 44 there is shown a data processing 
system 4400 upon which eiubodiments of the present invention, 
that is to say, the non-linear systems and the methods of 
design thereof can be implemented or realised. It will also 
be appreciated that the non-linear systems so designed may be 
implemented in digital form using the data processing system 
or other suitable hardware and software. The data processing 
system 4400 comprises a central processing unit 4402 for 
processing computer instructions for implementing the design 
of nonlinear systems given an input signal and particular 
output signal requirements. The data processing system also 
comprises a memory 4404 for storing data to be processed or 
the results of processing as well as computer program 
instructions for processing such data, a system bus 4406, an 
input device 4408, an output device 4410 and a mass storage 
device, for example, a hard disc drive 4412. 
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CLAIMS 

1. A method for designing a non-linear system for 
transferring energy from a time or spatial domain 
input signal having a first spectrum at a first pre- 
determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies, said method comprising the 
steps of 

identifying or specifying the first spectrum of the 
time or spatial domain input signal from which energy 
is to be transferred, 

specifying the second spectrum of the time or spatial 
domain output signal to which said energy is to be 
transferred, and 

calculating, using a frequency domain description of 
said output signal, for example, the output spectrum, 
expressed in terms of a frequency domain description 
of said input signal and coefficients of a time or 
spatial domain description of a generalised non- 
linear system, the coefficients of the time or 
spatial domain description of the generalised non- 
linear system in order to give effect to the energy 
transfer . 

2. A method as claimed in claim 1, further comprising the 
step of 

selecting a time or spatial domain description of the 
generalised non-linear systems- 
determining or defining a frequency domain description 
of the time or spatial domain input for the generalised 
non-linear system; and 
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determining or defining the frequency domain 
description of the output signal, for example, the 
output spectrum, of the generalised non-linear system 
expressed in terms of the frequency domain description 
of said input signal and the coefficients of the time 
or spatial domain description of a generalised non- 
linear system. 



3. A method as claimed in either of claims 1 or 2, wherein 
the frequency domain description of the input signal is 
U(jw), the time or spatial domain description of said 
generalised non-linear system is given by the 
generalised NARX model 

N 

y(k) = Zyn(k) (CD 



15 where 



p+q 



with 

p-hq = n, - 1, =1,...,K„, i = l,..„p + q, and J ^Z'-'Z 



(C2) 



the frequency domain description of the output of the 
generalised non-linear system is given by 

N 

Y(jw) = 2;Y„(jw) (C3) 

n=l 

where 

Yn(jw) = -i^^ K(jw,,...,jwjnU(jw,)do-, (C4) 

^^^^ w, + _, + w3=w i=l 

25 N is the maximum order of dominant system nonlinearities , 
J( Oda-, 
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denotes an integration over the nth-dimensional hyper- 
plane w^+,"v+w„ = w , and H^( jw^,---, jw^) , n = l,,,.,N are 
generalised frequency response functions of the non- 
linear system. 

4 . A method as claimed in any preceding claim, further 
comprising the step of 

determining a mapping between the time or spatial 
domain description of the generalised nonlinear system 
and the frequency domain description of the generalised 
nonlinear system. 

5, A method as claimed in claim 4, wherein the mapping from 
the time or spatial domain description of the 
generalised non-linear system to the frequency domain 
description of the system is given as 




= Z<=On(ll'---'ln) exp[- j(w,l,+,---,+W„l„)] 




q=l p=l li,lp,,=l 



)K-<,,p(jwi,---, 





(C5) 



where 



n-p+l 



•••,jw„) exp[- j ( 



-v+w,)lp] 



H„p(jw 



(C6) 
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6. A method as claimed in any preceding claim, further 
comprising the steps of 

defining or determining a general relationship between 
5 the input and output frequency or frequency ranges of 

the generalised non-linear system. 



10 



15 



20 



7 , A method as claimed in any preceding claim, wherein the 
relationship between the input and output frequencies 
or frequency ranges is given by the following 

fY=fY,Uf.,_, (C7) 
where f^ denotes the range of frequencies of the output, and 
fy_ and f^__^ denote the ranges of frequencies produced by the 
Nth-order and ( N -1 ) th-order nonlinearities, and 

nb 



when 
when 



(a + b) 

nb 
(a + b) 



na 
(a + b) 
na 



<1 



>1 



(C8) 



k=o vd-rjj; [_ (a + b) 

n = N and N-1 
where [,] relates to or means take the integer part, 

na 



1 = 



(a + b) 



+ 1 



25 



Ij^ ^[na-k(a + b) ,nb-k(a + b)] for k = 0,...,i* - 1, 
I. =[0,nb-r(a+b) ], 



and the frequencies of the signal to be processed are 
in the range [a,b] and given [a,b] and the required 
output frequency range fy^ the method further comprises 

the step of determining the smallest N from the 
relationship above for the generalised non-linear 
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system which can bring about the specified frequency 
domain energy transformation. 

A method as claimed in claim 7, wherein^ with N having 
been determined and K^,n = 1,...,N , being given a priori, 
the method further comprises the steps of: 

taking N as N and determining the coefficients of the 
time or spatial domain model of the generalised non- 
linear system from the description for the system 
output spectrum given in terms of the spectrum of the 
input signal and the coefficients of the time or 
spatial domain model of the said generalised non-linear 
system. 

A method as claimed in claim 8, further comprising the 
steps of 

substituting H^( j w^,-- jw^) given in (C5) into (C4), and 
substituting the resultant expression for Y^ijw) into 
(C3) to obtain the description for the system output 
spectrum in terms of a function of the spectrum of the 
input signal and the coefficients of the time or 
spatial domain model of the said generalised nonlinear 
system. 

0. A method for realising or manufacturing a non-linear 
system for transferring energy from a time or spatial 
domain input signal having a first spectrum at a first 
pre-determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies, the method comprising the steps 
of 
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(a) designing the non-linear system using the method as 
claimed in any of claims 1 to 9; and 

(b) materially producing the non-linear system so 
designed or using the non-linear system so designed to 
modify materially the transfer function of an existing 
linear or non-linear system. 

11. A data processing system for designing a non-linear 
system for transferring energy from a time or spatial 
domain input signal having a first spectrum at a first 
pre-determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies, said system comprising 

means for identifying or specifying the first spectrum 
of the time or spatial domain input signal from which 
energy is to be transferred, 

means for specifying the second spectriom of the time or 
spatial domain output signal to which said energy is to 
be transferred, and 

means for calculating, using a frequency domain 
description of said output signal, for example, the 
output spectrum, expressed in terms of a frequency 
domain description of said input and coefficients of a 
time or spatial domain description of a generalised 
non-linear system, the coefficients of the time or 
spatial domain description of said generalised non- 
linear system in order to give effect to the energy 
transfer . 

12. A system as claimed in claim 11, further comprising 
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means for selecting a time or spatial domain 
description of the generalised non-linear systems- 
means for determining or defining a frequency domain 
description of the time or spatial domain input for the 
generalised non-linear system; and 

means for deteimining or defining the frequency domain 
description of the output of the generalised non-linear 
system expressed in terms of the frequency domain 
description of said input signal and the coefficients 
of the time or spatial domain description of a 
generalised non-linear system. 

13. A system as claimed in either of claims 11 or 12, 
wherein the frequency domain description of the input 
signal is U(jw), the time or spatial domain description 
of the generalised non-linear system is given by the 
generalised NARX model 

y (k) = 2Iy,(k) (C9) 

where 

n K f^ p p + q 

yn(k) = s z Cp^(i,,...,ip,<^)n y (k - n u ( k - (CIO) 

P=01i,lp,q=l i«l i = p + l 

with 

p + q = n, li=l,...,K^, i = l,-..,p + q, and 21 = S*" S 



the frequency domain description of the output of the 
generalised non-linear system is given by 

N 

Y(jw) = SY„(jw) (Cll) 

n=l 

where 
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N is the maximum order of dominant system nonlinearities, 
J( .)da-„ 



5 denotes an integration over the nth-dimensional hyper-plane 
WjL+,--v+w^ = w , and H^( jw^,---, jw^) , n = l,...,N are generalised 
frequency response functions of the said non-linear system. 



14. A system as claimed in any of claims 11. to 13, further 
10 comprising 

means for determining a mapping between the time or 
spatial domain description of the generalised nonlinear 
system and the frequency domain description of the 
generalised nonlinear system. 



15 



20 



15. A system as claimed in claim 14, wherein the mapping 
from the time or spatial domain description of the 
generalised non-linear system to the frequency domain 
description of the system is given as 

|l- Z^iodi) exp[- j (w,-h,..v+w^) li]|H,(jWi,---,jwJ 

n-1 n-q Kp^^ 

+ Z E Z Cpq ( li ' • • ■ ' lp*q ) e xp - j(w „_^^, Ip,, +, • • • ,+w „ Ip,^ ) H„_c,p ( j w 1 , ■ • • , j w „_ J 

q=l p=l 1„V<,=1 *• J ^ / 



+ Z ZSo(ll'--'lpKp(^^l'"'"'^^n) 
p=2 1,,lp=l 

(C13) 

where 
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n-p+l 

(C14) 

t 

16. A system as claimed in any of claims 11 to 15, further 
5 comprising 

means for defining or determining a general 
relationship between the input and output frequency or 
frequency ranges of the generalised non-linear system. 



10 



15 



20 



17. A system as claimed in any of claims 11 to 16, wherein 
the relationship between the input and output 
frequencies or frequency ranges is given by the 
following 

fy^f^sUf.,.. (C15) 
where fy denotes the range of frequencies of the output, and 
fy. and fy_ ^ denote the ranges of frequencies produced by the 

Nth-order and ( N -1 ) th-order nonlinearities , and 

nb 



k=0 



when 
when 



l.k=0 



(a + b) 

nb 
(a + b) 



na 
(a + b) 
na 



<1 



>1 



(C16) 



(a + b) 
n=N and N-1 
where [.] relates to or means take the integer part, 

na 



1 = 



+ 1 



(a + b) 

1^ =[na-k(a + b) ,nb-k(a + b)] for k = 0,...,i' - 1 , 
I.. =[0,nb-i'(a + b) ], 



and the frequencies of the signal to be processed are in the 
25 range [a,b] and given [a,b] and the required output 



88 

SUBSTITUTE SHEET (RULE 26) 



wo 99/45644 PCT/GB9 9/00550 

frequency range fyr the system further comprises the means 

for determining the smallest N from the relationship above 
for the said generalised non-linear system which can bring 
about the specified frequency domain energy transformation. 



18. A system as claimed in claim 17^ wherein, with N having 
been determined and Kj^,n = 1,,.,, N , being given a priori, 
the system further comprises the means : 

for taking N as N and for determining the coefficients 
of the time or spatial domain model of the generalised 
non-linear system from the description for the system 
output spectrum given in terms of the spectrum of the 
input signal and the coefficients of the time or spatial 
domain model of the generalised non-linear system. 

19. A system as claimed in claim 18, further comprising 
means for substituting H^( j Wi,---, jw^) given in (C13) 
into (C12), and substituting the resultant expression 
for YJjw) into (Cll) to obtain the description for the 
system output spectrum in terms of a function of the 
spectrum of the input signal and the coefficients of 
the time or spatial domain model of the generalised 
nonlinear system. 



20. A computer program product for designing a non-linear 
system for transferring energy from a time or spatial 
domain input signal having a first spectrum at a first 
pre-determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies, the said product comprising, a 
computer readable storage medium comprising: 
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computer program code means for identifying or 
specifying the first spectrum of a time or spatial 
domain input signal from which energy is to be 
transferred, 

computer program code means for specifying the second 
spectrum of a time or spatial domain output signal to 
which said energy is to be transferred, and 

computer program code means for calculating, using a 
frequency domain description of the output signal, for 
example, the output spectrum, expressed in terms of a 
frequency domain description of the input and 
coefficients of a time or spatial domain description of 
a generalised non-linear system, the coefficients of a 
time or spatial domain description of said generalised 
non-linear system in order to give effect to the energy 
transfer . 

1. A computer program product as claimed in claim 20, 
further comprising 

computer program code means for selecting a time or 
spatial domain description of the generalised non- 
linear systems- 
computer program code means for determining or defining 
a frequency domain description of the time or spatial 
domain input for the generalised non-linear system; and 

computer program code means for determining or defining 
the frequency domain description of the output of the 
generalised non-linear system expressed in terms of the 
frequency domain description of said input signal and 
the coefficients of the time or spatial domain 
description of a generalised non-linear system. 
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22. A computer program product as claimed in either of 
claims 20 or 21, wherein the frequency domain 
description of the input signal is U(jw), the time or 
spatial domain description of the generalised non- 
linear system is given by the generalised NARX model 

N 

y(5<^) = Zyn(k) (C17) 

where 

p = 01i,lp,q=l i = l i = p + l 

with 

^ ^ ^ 

p-hq = n, 1^=1, ...,K„, i = l,...,p + q, and 2^ ^ X'" 2 



the frequency domain description of the output of the 
generalised non-linear system is given by 

N 

Y(jw) = 2Y„(jw) (C19) 

n=l 

15 where 

Yn(jw) = -i^ K(jw,,...,jwjnU(jw,)da-„ (C20) 

N is the maximum order of dominant system nonlinearities , 
J(.)dcr, 

20 denotes an integration over the nth-dimensional hyper-plane 
Wi+/-",+w^ = w , and Hj^( jw^^,---, jw^) , n = l,...,N are generalised 
frequency response functions of the non-linear system. 

23. A computer program product as claimed in any of claims 
25 20 to 22, further comprising 
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computer program code means for determining a mapping 
between the time or spatial domain description of the 
generalised nonlinear system and the frequency domain 
description of the generalised nonlinear system. 



24. A computer program product as claimed in claim 23, 
wherein the mapping from the time or spatial domain 
description of the generalised non-linear system to the 
frequency domain description of the system is given as 

|l- Z^iodi) exp[- j (w,+,.",+wJ lx]|H,(jw„...,jwJ 

= S^On(ll^---/ln)e^p[- j(Wili+,---,+W^lj] 

n-1 n-q K^^^ 

q=i p=lli,lp..i-l ^ /J ^ / 

p=2 1,Ap=l 

(C21) 

where 

n-p+l 

Hnp(jWi,---,jw^)= ZH,(jw,,--.,jwJH„.,^^,(jw,^,,---,jwJexpf-j (w,+,-..,+w,) 1 ] 

i=l L t'J 

(C22) 



25. A computer program product as claimed in any of claims 
2 0 to 24, further comprising 

computer program code means for defining or determining 
a general relationship between the input and output 
frequency or frequency ranges of the generalised non- 
linear system. 

26. A computer program product as claimed in any of claims 
20 to 25, wherein the relationship between the input 
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10 



and output frequencies or frequency ranges is given by 
the following 

f.= f.,Uf.,., (C23) 
where f^ denotes the range of frequencies of the output, and 
fy_ and fy__^ denotes the ranges of frequencies produced by 
the Nth-order and ( N -1 ) th-order nonlinearities , and 



k=0 



when 
when 



nb 



(a + b) 

nb 
(a + b)" 



na 



. (a + b) 
na 



<1 



>1 



(C24) 



k=0 V a -r ; L ( a + b ) 

n = N and N — 1 
where [ . ] relates to or means take the integer part, 

na 



1 = 



+ 1 



(a + b) 

I,, =[na-k(a+b) ,nb-k(a+b)] for k 0,,..,i* - 1, 
I.. =[0,nb-i*(a + b) ], 



and the frequencies of the signal to be processed are in the 
range [a,b], given ta,b] and the required output frequency 
15 range fy^ the computer program product further comprises 

computer program code means for determining the smallest N 
from the relationship above for the said generalised non- 
linear system which can bring about the specified frequency 
domain energy transformation. 

20 

27. A computer program product as claimed in claim 26, 
wherein, with N having been determined and K^,n = 1 N , 
being given a priori, said product further comprises 
computer program code means for taking N as N and 
25 determining the coefficients of the time or spatial 

domain model of the generalised non-linear system from 
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the description for the system output spectrum given in 
terms of the spectrum of the input signal and the 
coefficients of the time or spatial domain model of the 
said generalised non-linear system. 

5 

28. A computer , program product as claimed in any of claims 
2 0 to 27, further comprising computer program code 
means for substituting H^{ j w^,-- j w^) given in (C21) into 

(C20), and substituting the resultant expression for 
10 Y^(jw) into (C19) to obtain the description for the 

system output spectrum in terms of a function of the 
spectrum of the input signal and the coefficients of 
the time or spatial domain model of the generalised 
nonlinear system. 

15 

29. A method for designing a non-linear system for 
transferring energy from a time or spatial domain input 
signal having a first spectrum at a first pre- 
determinable frequency or range of frequencies to a 

20 time or spatial domain output signal having a second 

spectrum at a second pre-determinable frequency or 
range of frequencies substantially as described herein 
with reference to and/or as illustrated in the 
accompanying drawings . 

25 

30. A method for manufacturing a non-linear system 
comprising the steps of designing a non-linear system 
for transferring energy from a time or spatial domain 
input signal having a first spectr\m at a first pre- 

30 deteirminable frequency or range of frequencies to a 

time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies using a method substantially as 
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described herein with reference to and/or as 
illustrated in the accompanying drawings and materially 
producing the non-linear system so designed. 



31. A data processing system for designing a non-linear 
system 

for transferring energy from a time or spatial domain 
input signal having a first spectrum at a first pre- 
determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies substantially as described herein 
with reference to and/or as illustrated in the 
accompanying drawings . 

32 . A computer program product comprising computer program 
code means for designing and/or implementing a non- 
linear system for transferring energy from a time or 
spatial domain input signal having a first spectrum at 
a first pre-determinable frequency or range of 
frequencies to a time or spatial domain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies substantially as 
described herein with reference to and/or as 
illustrated in the accompanying drawings. 

33. A non-linear system for transferring energy from a time 
or spatial domain input signal having a first spectrxam 
at a first predeterminable frequency or range of 
frequencies to a time or spatial domain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies substantially as 
described herein with reference to and/or as 
illustrated in the accompanying drawings. 
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34. A non-linear system which can transfer energy from a 
time or spatial domain input signal having a first 
spectrum at a first pre-determinable frequency or range 
of frequencies to a time or spatial domain output 
signal having a second spectrum at a second pre- 
determinable frequency or range of frequencies, said 
system comprising 



or 



means for identifying the first spectrum of the time 
spatial domain input signal from which energy is to be 
transferred. 



or 



means for specifying the second spectrum of the time 
spatial domain output signal to which said energy is to 
be transferred, and 

means for giving effect to the energy transfer using 
coefficients of a time or spatial domain description of 
a generalised non-linear system, said coefficients 
having been calculated using a frequency domain 
description of said output signal, for example, the 
output spectrum, expressed in terms of a frequency 
domain description of said input signal and 
coefficients of a time or spatial domain description of 
a generalised non-linear system. 



A non-linear system for transferring energy from a time 
or spatial domain input signal having a first spectrum 
at a first pre-determinable frequency or range of 
frequencies to a time or spatial domain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies, said system 
comprising means for implementing a method as claimed 
in any of claims 1 to 9 or 29, 
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36. A non-linear system for transferring energy from a time 
or spatial domain input signal having a first spectrum 
at a first pre-determinable frequency or range of 
frequencies to a time or spatial domain output signal 
having a second spectrum at a second pre-determinable 
frequency or range of frequencies, said non-linear 
system comprising a data processing system as claimed 
in any of claims 11 to 19 or 31. 

37. An article of manufacture comprising a computer usable 
medium with computer readable program code means 
embodied in the medium for designing a non-linear 
system for transferring energy from a time or spatial 
domain input signal having a first spectrum at a first 
pre-determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrim at a second pre-determinable frequency or 
range of frequencies, the computer readable program 
code means in said article comprising: 

computer readable program code means for identifying 
the first spectrum of a time or spatial domain input 
signal from which energy is to be transferred, 

computer readable program code means for specifying the 
second spectrum of a time or spatial domain output 
signal to which said energy is to be transferred, and 

computer readable program code means for calculating, 
using a frequency domain description of the output 
signal, for example, the output spectrum, expressed in 
terms of a frequency domain description of the input 
and coefficients of a time or spatial domain 
description of a generalised non-linear system, the 
coefficients of a time or spatial domain description of 
said generalised non-linear system in order to give 
effect to the energy transfer. 
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(a) 




Given a sampled series {u(k)} of the signal | 
u(t) tp process and the sampling interval T | 




r 


Given a range [c,d] of frequencies specified 
for the output components of the filter to design 




r 


Given a spectrum Y*(jw) s 
of the filter output over t 


Specified for the spectrum | 
he frequency band [c,d] | 



3002 



3004 



3006 



(b) 



Calculate U (jw), the spectrum of the input, using 
a fast Fourier Transformation algorithm to yield 

Hf-) ° f 

where M is the length of the sampled series {u(k)} 



Evaluate the range [a,b] of the frequencies in u(t) from 
the computed input spectrum 



Setn=1.No=0 I 3012 




Perform the following computations 
na 



(a+b) 



+ 1 



Ij=[na-i(a+b). nb-i(a+b)] 
for i=0,....i*-1 
ljv=[0,nb-i*(a+b)] 



0 

Fig. 30 
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© 



Evaluate the frequency range f w produced 

^n 

by the nth-order nonlinear output as 
"Uli when ' ™ 



i=1 
i' 



[J Ij when 



U=0 



(a+b) 
na 

(a+b) 



(a+b) 

na 
(a+b) 



<1 



>1 



where [ x ] denotes the integer part of x. 



(^^Whether A/q= 0 ?^ 3018 




3016 



Whether [c,d] has completely fallen into i^^^y '-'V 



3024 



Yes 



N=n I 3027 



SetK^=1. forn=NQ N 



© 



3028 



n 



Calculate 

f exp [-j(w^ l^+ n U (jw^) da^ 

J i=1 
w^ + +w^=w 

numerically at w=27t i/MT, i=0 M/2. 

for 1^ =1 1^=1 K^, and n=NQ N 




3030 
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(d) 



(e) 



L 3044 




Compute 
L= round 



cMT 

27C 



and i . = round 
d 



dMT 
2n 



where round (x) means to take the 
integer nearest to x 



3032 



K 



Based on equations 

n=No(2^) 1^.1^=1 n 

f exp[-j(w^1^+ +Wn1n)n Uaw^)dc^ 

•' ... PIT . i=1 



'='c-'c+'--'d 



Compute 

^On^'r - ' 'n)' 'l=' Kn....,in=''- -'^n- n=No, -.N. 

using at least square routine. 



3034 



Design a conventional linear filter H(q'^) to improve the 
above design result within the frequency band [c.d] and 
attenuate the responses over the frequency bands which 
are beyond [c.d]. 



3036 



set Kn =Kn+1. 
for n=N^ N 



(f) 



Check the frequency domain r^ponse 
of the filter ^ K, 

y(k)=H(q-1)=^ ^ Cq^Oi In).nu(k-I) 

over the frequency band [c,d] 



3038 



No 



Whether the response, in comparision 
with the specified response Y* (jw), 

is okay? ^.i-.-^040 

Yes 



The design is finished and the discrete 
version of the designed nonlinear system is 

N ^' 



n 



y(k)=H(q-1 ) ^ co^di In). U u(k-lj) 




3042 



End 
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(a) 



(b) 




Given a sampled series {u(k)} of the signal 
u(t) tp process and the sampling intenaal T. 



Given a range [c,d] of output frequencies 



H* 

Given a bound Y° (w) specified for the spectrum 

of the filter output over the frequency band [c,d] 



Calculate U (jw). the spectrum of the input, using 
a fast Fourier Transformation algorithm to yield 



where M is the length of the sampled series {u(k)} 



0 M 

2 



Evaluate the range [a.b] of the frequencies in u(t) from 
the computed input spectrum 



Setn=1, Nq=0 I 



3112 




Perform the following computations 
na 



+1 



(a+b) 

I j =[na-i(a+b). nb-l(a+b)] 

fori=0 i*-1 

l..=[0.nb-i*(a+b)] 



0 

Fig. 31 
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1 



Evaluate the frequency range fy produced 

^n 

by the nth-order nonlinear output as 
U Ij when 



'n 



i=1 



(a+b) [ (a+b) 



na 



1) Ij when 
i=0 (^■^'^^ 



na 



<1 



>1 



(a+b)_ 

where [ x ] denotes the integer part of x. 




3116 



Whether [c.d] has connpletely fallen into t,=fv Ufw 



3124 



Yes 



N 



=n] 3127 



Calculate |U|*...*|U (jw)| , the nth-order convolution 



n 



integration on |U (jw)| , numerically to yield 
|U|*-*|U (j27l i/MT)| , i=0 M/2. for n= A/q ,...,N. 



3128 
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Compute 






'dMT" 




round 


'cMT' 






and =round 


. 27C , 












where round (x) means to take the 




integer nearest to x 





3132 



Based on equations 

Yli^^i^)^n iU|-...-|u[i(f,i) 
n=NQ ^'^ • —i. 



n 



Compute 

Cp, n=NQ N, 

using at least squares routine under the constraint that 
the results must be positive 



3134 



Select NARX model parametres 

•^n' ^on^'l 'n) • 'l ^n'-"^ 'n=^ ^n* 

for n=NQ....,N, under the constraints 



Kn 



4=1 1^=1 



n=NQ N 



and other constraints specified by specific 
design problems 



3136 



Checl^ the frequency domain response 
of the filter ^ 

y(k)=H(q-1)=]£ 2^ Condi Ip)..!! 



u(k-L) 



n=Nn L I =1 
over the frequency bands beyond [c,d], 
which should, ideally, be negligible 



3138 




Whether the frequency response 
beyond the band [c,d] is ok? 



3140 
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© 



(f) 



The design is finished and the discrete 
version of the designed nonlinear system is 

N "^n n 
y(k)=H(q-1)=Y T Condi In)..llu(k-ip 



3142 




Design a conventional linear band pass filter 
H ( q-l) to attenuate the responses over the 
frequency bands which are beyond [c,d] 



3144 



Complete the design by constructing a nonlinear 
filter the discrete version of which is 

K, 



'n 



y(k)=H(q-1)=2^ 2^ Cq^Oi In)..riu(k-L) 




3146 
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